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Abstract 



The vacuum correlations of the gravitational field are highly non-trivial to be defined 
and computed, as soon as their arguments and indices do not belong to a background but 
become dynamical quantities. Their knowledge is essential however in order to understand 
some physical properties of quantum gravity, like virtual excitations and the possibility 
of a continuum limit for lattice theory. In this review the most recent perturbative and 
non-perturbative advances in this field are presented. (To appear on Riv. Nuovo Cim.) 
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1. INTRODUCTION. 



The aim of this paper is that of giving an introduction to the issue of vacuum corre- 
lations in quantum gravity and a review of the main advances recently made in this field. 
The whole subject is quite new and in rapid development, so the present account can only 
be considered as provisional. 

From the physical point of view, this is a very interesting field of research, since it 
concerns the physical modes of propagation of the gravitational field and their quantization, 
the structure of the vacuum state and the general problem of the observables in quantum 
gravity. 

It is known that a completely consistent theory of quantum gravity has not been 
established yet. Some of the autors we shall mention in this paper really work in this 
direction: the simplicial quantum gravity of Hamber (Chapter 7), for instance, aims to put 
Einstein's gravity on a solid quantum basis by simulating numerically its non-perturbative 
behaviour; also the original purpose of the work of Tsamis and Woodard (Chapter 2) was 
that of finding the "true" observables of quantum gravity, free of ultraviolet divergences. 

Our philosophy throughout this paper, however, will be mainly that of regarding 
quantum gravity as an effective quantum field theory, which has General Relativity as its 
classical limit, but which could go over to some more fundamental theory at very short 
distances. Our approach will be in some sense "phenomenological" , as we believe that 
knowledge of the properties of the gravitational vacuum correlations and other observables 
can help in guiding the fundamental research. 

From the technical point of view, the issue is confronted either in a suitable physical 
gauge (the Mandelstam covariant and the radial gauge, Ch.s 2, 3, 4) and in gauge-invariant 
form (Ch.s 5, 6, 7). In the following of this Introduction, after giving some general defi- 
nitions, we shall survey both approaches giving an overall outline of the paper. A more 
detailed outline for each chapter can be found at the beginning of the chapter itself. 

1.1 Vacuum correlations in the absence of gravity. 

The vacuum correlation functions are very important quantities in quantum field 
theory. It was shown by Wightman many years ago [Wightman, 1956] that the knowledge 
of all the correlation functions of a scalar field 4>{x) gives a complete information about 
the quantum dynamics of the field. More recently, in the lattice formulations of euclidean 
quantum field theory, the field correlations characterize the "phases" of the equivalent 
statistical system. Typically, they behave like (<j)(x)<f)(y)} ~ e -\ x ~v\/^ outside the region 
of phase transition, and like ((f)(x)(f)(y)) ~ \x — y\~ m in correspondence of the transition 
(where the continuum limit is recovered; see Chapter 7). 

In general, for vacuum correlations of a quantum field theory whose fields are denoted 
by $>a(x), &b(x) ••• we mean the expectation value on the vacuum state of the product 
of the fields: 




(1.1) 
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In this expression the indices A, B, ... are Lorentz indices (possible "internal" indices are 
understood) and the coordinates the coordinates of Minkowski space. 

In the case of free fields - for instance, scalar fields - the function G(x, y) is computed 
exactly and is called the Pauli- Jordan distribution 

(0(*)0(</)>o = ^53 / d 4 ke-^e(k°)5(k 2 -m 2 ). 

From this expression one can verify that also in a "trivial" system like a free field the 
vacuum correlations exibit interesting properties: for instance, they do not vanish when 
(x — y) is spacelike, although all the observables in x and y commute in this case. 

Intuitively, from the statistical point of view, the product of two real random variables 
vanishes on the average if the two variables are independent, and have thus the same prob- 
ability to assume sign + or -, with various amplitudes. But since the field $ propagates, 
its value in y depends on that in x, and averaging on all the configurations like in (1.1) we 
find that the correlation increases like \x — y\~ m when the two points get closer and closer. 
Of course, in the functional integral (1.1) all the field configurations are contained, not 
just the classical configurations which minimize the action; but if $ oscillates too rapidly 
at small distances in some configuration, the kinetic term of the action is very big and that 
configuration is strongly suppressed. 

As we saw before, the function G possesses in general, in the absence of gravity, some 
Lorentz indices. For instance, for the correlation of two vector fields we have 

Gp V (x,y) = (V^(x)K(y))o- 

This means that the field V^x) and the field V l/ (y) can be "rotated" by a Lorentz ma- 
trix - the same in x and in y - and the correlation function also transforms correspond- 
ingly; for instance, the correlation (Vi(x)V2(y))o becomes, in a rotated reference frame, 
(V'2(x')V's(y'))o- What allows us to compare in a consistent way two components of the 
field in two distinct points is the "rigid" structure of Minkowski space. Also it is obvious 
- still in the absence of gravity - that the distance \x — y\ is a fixed quantity, that does 
not depend on the field V itself. 

If we allow a local invar iance in our theory, the situation changes remarkably. This can 
be seen already in the usual gauge theories, based on "internal" symmetries. In this case, 
the two fields A^(x) and A u (y) which have to be correlated gauge-transform independently 
in x and y. This reflects itself in the fact that the correlation function is gauge-dependent. 
However, the coordinates x and y and the Lorentz indices are still referred to a rigid 
background. 

In the case of gravity, the analogue of the local gauge transformations are the coordi- 
nates transformations. There is no fixed background and the only way to give a physical 
meaning to a correlation function (except for the pure topological content) is that of com- 
paring the fields in two different points by parallel transport. As the distance of the points, 
the geodesic distance must be considered. But either in the evaluation of the parallel trans- 
port and of the geodesic distance the same dynamical field enters, whose correlations we 
are looking for. It is clear that this will introduce in general new technical difficulties. 
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1.2 Gravitational correlations in a physical gauge. 

If we are interested in all the "components" of the correlation functions, we have to 
compute them in a suitable preferred physical gauge. A first attempt to the computation of 
this kind of "physical Green functions" for gravity was done by Tsamis and Woodard (Sec- 
tion 2.2). They computed the vacuum correlations of the Mandelstam covariant (Section 

2.1) to one-loop perturbative order, in the hope that these quantities, being in principle 
observable, would have been free of divergences. Unfortunately, as a result of their hard 
work (one of the biggest perturbative computations in quantum gravity) they found that 
the usual divergences are still there, and some new ones appear, which are peculiar of the 
quantity under consideration. 

On the other hand, the Mandelstam formalism cannot be considered in general as a 
gauge-fixing and thus many standard techniques of QFT cannot be applied to this case. 
Independently of Tsamis and Woodard, Toller [1988] suggested to try to define the fields 
themselves with reference to geodesies and parallel transport, in order to give a direct 
physical significance to their correlation. This led to the classical radial gauge (Chapter 
3), which, on account of its simplicity and of the field-potential inversion formula (Section 

3.2) , has proved to be useful for other purposes too (Sections 3.4, 3.5). 

In this gauge the fields are defined in such a way that the vierbein e^(£) in a point 
£ coincides with that which is parallel transported from the origin to £ along a straight 
line. Moreover, this line is the (locally unique) geodesic which connects the origin to the 
point £ and the geodesic distance from the origin is given by |£|. Thus in this gauge the 
correlation function (e^(0)e^(£)) automatically satisfies the criteria we enunciated at the 
end of Section 1.1. 

Inserting the radial gauge fixing in some version of the functional integral of gravity, 
one may compute the corresponding propagator. This propagator gives automatically the 
lowest-order physical correlation function and allows in principle to compute this correla- 
tion to any order. 

As expected, the propagator in the radial gauge is not trivial (including in fact the 
geodesies and parallel transport) and contains some singularities (Sections 4.2, 4.3). By 
means of a suitable projection procedure, however, it is possibile to gain an insight into the 
singularities and to eliminate them. It is believed now that the "additional" divergences 
which appear in the mentioned work of Tsamis and Woodard could be eliminated through 
a similar procedure. There is no special connection between these divergences and the 
usual ultraviolet divergences of quantum gravity. 

Summarizing, we have in the radial gauge a complete solution to the problem of 
computing perturbatively the vacuum correlations at geodesic distance. The limitations 
are those inherent to perturbation theory in quantum gravity (i.e., non renormalizability 
for the Einstein theory, or non unitarity for the (R + _R 2 )-theory), and those due to the 
mathematical complexity of the radial propagator. On the other hand, some problems 
typical of algebraic gauges, like the need for special prescriptions in order to define uniquely 
the propagator [see for instance Bassetto, Nardelli and Soldati, 1991; Gaigg and Kummer, 
1990], do not affect the radial gauge. Also, in four dimensions the four-particles interaction 
vertex is likely to vanish in radial gauge since the four fields should be at the same time 
orthogonal to £ and between themselves. 
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1.3 Invariant gravitational correlations. 

The formal advances in the determination of the radial propagators do not say us 
much on some physical features which characterize the gravitational vacuum correlations. 
In particular, we may ask: which are the main differences with respect to the other gauge 
fields? Which are the modes of the field which are really correlated? Is the gravitational 
interaction truly due to the exchange of virtual particles, and thus connected to the vacuum 
fluctuations in the same way like the other fundamental interactions do? 

In Chapters 5 and 6, the first attempts to answer these questions are presented. 
The crucial quantities under investigation are the Wilson loops of the connection, which 
are scalar quantities and give a sensible information about the geometry of spacetime. 
The computations we present here start from a flat background and take into account 
corrections up to order %. This is justified on length scales much larger than Planck's 
length; at smaller distances, non-perturbative methods must be employed (see Chapter 7) . 

In Chapter 5 the Wilson loops of the gravitational connection are defined (Section 
5.1) and computed to order % in Einstein's theory (Section 5.4). The invariant two-point 
functions of the Riemann curvature are computed too (Section 5.3). 

The Wilson loops vanish and this peculiar result, when analyzed in terms of the rele- 
vant gauge group for the euclidean theory, 50(4) (Section 6.1), shows that the functional 
integral of quantum gravity does not contain to order h any configuration with localized 
curvature. Such a behaviour contrasts strongly with that of usual gauge fields, and may 
be interpreted as the absence, to this order, of true virtual gravitons (those which appear 
in the diagrams for the cross sections or similar represent excitations of the metric which 
do not carry curvature!). It can be shown that this result holds more generally than just 
for Einstein gravity, but relies on the flat background. 

Chapter 6 contains also some more material on the Wilson loops: the correlation 
function of two loops to leading order and its "geodesic corrections" (Sections 6.2, 6.3). In 
Section 6.4 the formula for the static potential energy in quantum gravity is presented. 

Finally, Chapter 7 is dedicated to lattice gravity, with special regard to the "quantum 
Regge calculus" of Hamber. The outcome of the recent Montecarlo simulations is very 
promising: the transition between the two phases of the lattice system appears to be of 
second order, and this has been confirmed by the behavior of the (simple) correlation 
functions which have been computed up to now. 

Notations. 

We refer mainly to spacetime dimension 4; otherwise the dimension is denoted by N 
and assumed to be bigger than 2. The units are such that c = h = 1; in some formulas 
where the classical contributions need to be distinguished from the quantum ones, h is 
written explicitly. The symbol k (or sometimes £) denotes Planck's length \J1Q-kG. In some 
of the reviewed papers the metric has euclidean signature, in other it is minkowskian; we 
have used both notations, although for gravity the correspondence between the euclidean 
theory and the theory on Minkowski space is not completely established yet [see for instance 
Mazur and Mottola, 1990]. 
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The derivatives of the form are taken with respect to the whole argument of the 
function. A suffix like [ab] denotes antisymmetrization (without the factor |). The notation 
"o( )" means "of higher order than"; "0( )" means "of the same order as". 

For recent general reviews on quantum gravity and references see for instance Alvarez 
[1989] or the lectures of the Les Houches Summer School, 1992. 
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2. THE MANDELSTAM COVARIANTS. 



In this first, "historical" chapter we shall first introduce the general classical Man- 
delstam covariants according to the original definition [Mandelstam, 1962], with a special 
emphasis on their geometrical meaning (Section 2.1); then we shall present in short the 
perturbative computation, due to Tsamis and Woodard, of a special kind of two-point 
functions inspired to Mandelstam's covariants (Section 2.2). 

2.1 The original definition. 

In the early Sixties it was already clear to a few quantum theorists that in a true 
theory of quantum gravity the gravitational field had to play the hard role of furnishing 
coordinates to spacetime while being at the same time a quantum object. 

In the mentioned paper, Mandelstam dicussed a technique for defining spacetime 
and its point without using coordinates, but replacing them with "paths" , which may be 
themselves influenced by the (eventually quantized) gravitational field. 

His pioneering work, although technically heavy in some points, had a great influence 
on the subsequent developments of gauge theories. Also, it is very remarkable that through 
his paths-based formalism Mandelstam succeeded in [1968] in finding the correct Feynman 
rules for gravity and Yang-Mills theories; these rules were derived at the same time by 
De Witt [1967 b, c] and by Faddeev and Popov, whose functional-integral technique has 
become later the standard one. 

In the Mandelstam formalism the coordinates of spacetime do not appear. A point 
is defined as the end of a path, which in turn is specified by an infinite list of "small 
steps" to be taken in sequence by an observer, with reference to his local frame. After each 
infinitesimal step, the frame is parallel transported to the new location, and so on. 

Thus the fields are not functions of four coordinates (xq, x\, X2, £3), but functionals 
of "paths" P. On the other hand, the paths P are purely geometric, intrinsic objects, 
while the coordinates (xo, xi, x%, X3) are subject in General Relativity to a vast class of 
tranformations without change of the physical results. 

The basic gravitational field is the anholonomic Riemann curvature R a bcd(P)- The 
components of R are referred to the parallel transported local reference frame, and are 
then coordinate-independent. The fields can be canonically quantized and their correlation 
functions defined, through a pretty complicated procedure. 

In a recent work of Teitelboim [1993] a more accurate mathematical definition of 
the classical Mandelstam's theory is given, also showing how the field equations in path 
space may be derived from an action principle and how is it possible to decide consistently 
"when two paths end at the same point" - a concept which was undefined in the original 
Mandelstam's theory. 

An interesting possibility within the Mandelstam formalism is that of fixing the paths 
through some prescription. A natural prescription, corresponding to straight lines in flat 
spacetime, is that of taking geodesic paths. In fact, the direction of a geodesic line, as 
observed in the free- falling local system, is constant, so the "infinite" list defining the 
path P is very much simplified. This geometrical construction is explicitly implemented 
- starting from a given central point - in the radial gauge (Chapters 3, 4), which thus 
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establishes a connection between a special case of Mandelstam formalism and the usual 
formalism of General Relativity or Einstein-Cartan theory. 

2.2 The computation of Tsamis and Woodard. 

The starting point of this work [1992] was that of considering the correlation of two 
Mandelstam-like fields, placed at the ends of a geodesic line of fixed length. This corre- 
lation was written by inserting the suitable matrices of parallel transport and evaluated 
perturbatively to one-loop order. 

The motivation for this computation originated in the troubles of quantum gravity 
with the unrenormalizable ultraviolet divergences. Some parallels with QED suggested that 
these divergences could cancel when the true "physical Green's functions" are computed, 
instead of generic quantities which are subject to coordinates transformations and are thus 
unphysical. 

The paper of Tsamis and Woodard is very remarkable from the "phylosophical" point 
of view and for its quantum- field content. 

"Instead of changing the fundamental dynamical principle", they proposed "to alter 
the way in which physical questions are asked. This idea is hardly new, nor is it unique 
to quantum gravity [...]. The particular detail we propose to treat more carefully is the 
physical coordinate system. We shall do this by introducing a new interpolating field, the 
Mandelstam covariant [...]. The Mandelstam covariant is not a local, invertible redefinition 
of the usual field [...]. It follows that the Green's functions deriving from the Mandelstam 
covariant need not interpolate the usual S-matrix, nor contain the usual ultraviolet diver- 
gences. Indeed we will prove that the usual result does change. Unfortunately it does so 
in the direction of greater divergence; however, the reason for this seems to be correctable 



And going on with the citation: "The issue for us is not so much that local Green's 
functions depend upon the gauge but rather that this gauge dependence makes them 
unwieldy experiments. In fact none of them is a reliable probe of quantum gravity, because 
they call for measurements that can never actually be made [...] and must we be very 
concerned over the fact that the theory predicts a divergent result for them?" 

Tsamis and Woodard write the Mandelstam covariant as 



K abcd [e](V, x) = M k a [e](V : x)M l b [e}(V 7 x)M?[e](V, x)M%[e](V, x)R klmn ( X [e](l, V", x)) 



where the bar over the Riemann tensor denotes that its indices are local Lorentz: 



By x[e](r, V, x) the geodesic path is meant which carries the Mandelstam observer from 
his origin x to the point at which he measures the curvature, r is the geodesic parameter 
and V the initial direction. It is always assumed that the fundamental dynamical variable 
is the vierbein e. M is the matrix of the parallel transport: 




eaa(z)4(z)e2(z)e 8 d (z)R^ s (z). 



M b a le](V,x) = e a ce (x(l))M$le](V,x)e?(x) 
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It can be proved that 1Z can not be related to R by any local field transformation. 
The reason for this is essentially the appearence in 1Z of the non-local matrix M. 

The one-loop computation of the correlations of TZ is extremely long. The expectation 
values are organized into 30 "/c-point functions" which were evaluated using the dimensional 
regularization. The "^/-ordering" corrections (corresponding to the path ordering) were 
then added, and the final result expressed as a linear combination of 8-index "master tensor 
structures" . 

At each order the result breaks up into "uncorrected terms" (corresponding to the 
naive flat space perturbations), "length corrections", "index corrections" and "length and 
index corrections". This is precisely what is expected from the physical features of the 
gravitational correlation functions (see the Introduction, and Section 6.3 for a concrete 
example) . 

The "additional singularities" which affect the result have the same form like those 
of the "(P° P°) propagator" of radial gauge. For this reason, it is very likely that they 
could be eliminated, in principle, by a suitable projection procedure (compare Section 4.2, 
4.3). In practice, however, the algebraic structure of the Mandelstam covariant employed 
by Tsamis and Woodard is too complicated to allow this. For a comparison with the 
radial-gauge formalism, see the beginning of Section 4.4. 
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3. CLASSICAL RADIAL GAUGE AND APPLICATIONS. 



In this chapter an introduction to the classical radial gauge is given. In Section 3.1 
the gauge condition is defined (in the vierbein formalism) and its attainability is shown. 
In Section 3.2 we write the formulas which allow to express the radial vierbein and the 
radial connection in terms of the Riemann tensor and of the torsion (the so-called inversion 
formulas). In Section 3.3 the transformation properties of the radial fields under rotation 
or translation of the central vierbein are given. In Section 3.4 we illustrate in short a 
remarkable application of the radial gauge to (2+l)-gravity and finally, in Section 3.5, 
an application to the description of the geodesic motion of test particles in a fluctuating 
gravitational field. 



3.1 Main features. 



The radial gauge x^A^(x) = has been proposed many years ago by Fock [1937] and 
Schwinger [1952, 1989; for a complete bibliography see M. and Toller, 1990] for electrody- 
namics or Yang-Mills theory. Its most interesting property is the so called fields-potentials 
inversion formula. Let us suppose that the field strength F fJiU (x) is known. In order to 
compute from F^ u (x) the vector potential A^(x) in radial gauge, one must solve the system 

x^A^x) = 

d^A v {x) - d v A„{x) + ig[Ap{x),A v {x)] = F^{x). 
Its general solution is 

A ft (x) = C d\\x v F^{\x) + A* om {x), 
Jo 

where the homogeneous part A^ om (x) is solution of the system 
x^ om (x) = 

d»A* om {x) - d v A* om {x) + ig[AX° m (x),AX° m (x)} = 0; 

that is, we have 

A "° m[l) = ^ w 

being / an arbitrary homogeneous function of degree 0. If one wants A^(x) to be regular 
at the origin, / must be set equal to zero. In this case, the gauge condition can be regarded 
as "complete" , in the sense that no residual gauge is present. Until the introduction of the 
general radial projectors (Section 4.2), we shall stay here in the hypotesis that all fields are 
regular at the origin; thus our inversion formula will be simply given by the "0-1 projector" 

A^x)= [ d\\x v F^(\x). (3.2) 
Jo 
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Let us now consider the case of gravitation. We recall that in order to describe the 
gravitational field it is sometimes convenient to use, instead of the metric tensor g pv {x), 
the vierbein fields e p (x) and e%(x), which are defined by 

e^(x)el(x)S a b = 9nu(x); (3.3) 
e£(x)el(x) = 6*; eS(x)e»(x) = S b a . (3.4) 

The holonomic indices n, u,... and the "internal" indices a, 6, ... range, in general, between 
1 and N, where N is the dimension of spacetime; the N vectors {ei(x), ejv(a;)} can be 
thought to represent a local reference frame at any point x. Condition (3.3) determines 
e p (x) up to a local rotation; the gauge potentials T a p and field strengths R ab v corresponding 
to this local invariance are related to the Christoffel symbol and to the Riemann tensor by 
the formulasf 



rZ, = rt a e b p e£ + eZd a e a p] (3.5) 
K pa = e»elRt pa . (3.6) 

In this formalism, two gauge fixing conditions are needed, in order to eliminate the 
freedom to make diffeomorphisms (i.e., coordinate transformations), and the freedom for 
local rotations of the vierbein. Recently a radial gauge condition has been proposed also 
for gravity. It has the form [M. and Toller, 1990] 

£T^(£) = 0, (3.7a) 
£"e»(0 = £"^- (3-76) 

The gauge-fixing conditions (3.7a) and (3.7b) have a simple geometrical interpretation. 
Condition (3.7a) means that the tetrads are parallel transported from the origin along the 
straight lines of the form 

R, < s < 1}; (3.8) 

condition (3.7b) means that these lines are auto-parallel, i.e. they are geodesic lines (in 
the absence of torsion). The coordinates £ coincide then with the well-known "normal 
coordinates" [see for ex. Kobayashi and Nomizu, 1969]. 

It is possible to give explicit formulae for the calculation of the normal coordinates 
and of the parallel-transported tetrads of the radial gauge. To this end, we consider an 

f The connection V? is viewed as a variable independent of the vierbein; so this for- 
malism is often referred to as "first order", while the metric formalism is defined to be 
of the "second order" . As we shall illustrate below, it is possible to define a radial gauge 
condition also in the second order formalism. We do not dwell here upon the relationship 
between the two formalisms and the role of the torsion [see, for instance, Hehl et al., 1976]. 
We notice, however, that the radial gauge in the first order formalism is meaningful also in 
spaces without metric connection and for gauge groups different from the Poincare group. 
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arbitrary coordinate system x^ and we indicate by Xq the coordinates of the origin P . We 
assume that the holonomic coefficients of the connection T^ p (x) are known. We indicate by 
the components of the tetrads with respect to the natural holonomic basis determined 
by the coordinates x^. The components of the dual tetrads in the natural basis determined 
by the normal coordinates £ p can be computed by means of the formula 

Since the tetrads are parallel-transported along the line (3.8), we have 

dt b {\) 



dX 



-KAx(X)]etWe T a (X)v a (3.9) 



and the fact the the vector with constant anholonomic components v a is tangent to the 
line (3.8) is expressed by the equation 

*^ = ett\)v*. (3.10) 

The equations (3.9) and (3.10), with the initial conditions :r M (0) = Xq and e^(0) = 
e^(-Po); determine the quantities x^ and as functions of A and of v a . However, it is 
easy to see that they depend on a particular combination of these variables, namely on the 
normal coordinates ^ a =Xv a . If we are considering a metric space and a metric connection, 
we can choose an orthonormal initial tetrad e o (0) and all the parallel-transported tetrads 
are automatically orthonormal. 

It is convenient to transform these differential equations with their initial conditions 
into the following pair of coupled integral equations 

^(0 = < + T f'e^OdX, (3.11) 
Jo 

= <(0) - T f K T [x(Xi)] e£(A£) e^(A£) dX. (3.12) 
Jo 

It is possible to solve these equations perturbatively to any desired order in f, by substi- 
tuting at each step the lower order solution in the left hand side integrals. This shows that 
the gauge condition is "attainable" in the sense that given the fields in an arbitrary gauge, 
it is always possible to compute the corresponding fields in the radial gauge. A simple 
application of these equations is the computation to lowest order in T of the coordinate 
transformation :r M (£). One obtains 

* M (£) = e - er f dt (i - 1) t»M) + o(f 2 ). (3.i3) 

Jo 

This formula will be used in Section 3.5 and 6.3. 
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Note that if we change the initial conditions by performing a Lorentz transformation 
L of the tetrad e a (0), (or more generally a transformation of the gauge group) we get a new 
solution obtained from the old one by means of the same Lorentz transformation L applied 
to all the tetrads and to the normal coordinates. This property is not trivial, because the 
parallel transport is generally non-commutative. A translation of the origin P affects the 
normal coordinates and the parallel-transported tetrads in a more complicated way, which 
will be illustrated in Section 3.3. 

The radial gauge conditions (3.7a) and (3.7b) can be regarded, in a sense, as an op- 
erational prescription which permits one to locate the measuring instruments in a neigh- 
bourhood of the observer, who lies at the origin P . In fact a simple way to explore this 
neigbourhood is to send from the origin many "space-probes" carrying clocks, gyroscopes 
and measuring instruments. A space-probe will be launched with four-velocity v a with re- 
spect to the given tetrad e a (0) and, if r is the proper time measured by the clock, £ a = rv a 
are the normal coordinates (in the absence of torsion). Of course, in Minkowski spacetime 
only the interior of the future cone can be explored in this way. 

Any space-probe will be able, in principle, to measure by local experiments the an- 
holonomic components of the fields at any point £. Moreover, let us go on with our 
Gedankenexperiment and suppose that any space-probe emits all the time some peculiar 
light signals along the axis of its local frame; the signals are received and recorded by the 
neighboring space-probes. By collecting all the data, the mentioned "central observer" will 
thus be able to compute also the vierbein e®(£). 

Finally, we remind that the radial gauge can be introduced also in the so-called second 
order formalism, that condition on the metric tensor [Menotti, M. and Seminara, 

1993]. In this case it is defined by 



A number of properties of the radial gauge in the tetrad formalism still hold in the 
metric formalism. We enumerate them without proof: 

- if <7 MI ,(£) is regular at the origin, then we have g^v(0) = r]^ u ; 

- coordinates £ satisfying (3.14) are normal coordinates; 

- the gauge is attainable, under some regularity hypotesis. 

An equation analogous to the inversion formulas (3.17), (3.19) holds just in the lin- 
earized theory. If we denote by R L the linearized Riemann tensor, we have 



(3.14) 




provided \R L 



(x)\ < \x 



-2+e 



dbS X 



o. 



3.2 Inversion formulas. 
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We shall now derive from the radial gauge conditions (3.7a) and (3.7b) two formulas 
analogous to (3.2), which give the field potentials r£ (£) and e£(£) in terms of the Riemann 
tensor RfnviO an< ^ the torsion tensor S^ v (^). 

We would like first to make the following observation. Remember that in this chapter 
we work in the hypotesis that the fields are regular (and differentiable) at the origin. Thus, 
taking the derivative of eq.s (3.7a) and (3.7b) we obtain 

r^(o) = o, 4(0) = *}. 

These equations show that in the radial gauge the gauge potentials at the origin take the 
values they have in a flat space. In other words, it is possible to eliminate the gravitational 
field at a given point. This is a formulation of the equivalence principle which is valid also 
in the presence of torsion, when it is not possible to eliminate the holonomic connection 
coefficients at a given point [see the references in M. and Toller, 1990]. 

Actually, one of the motivations which led to the formulation of the radial gauge was 
the need of a generalization of the equivalence principle to 10-dimensional theories defined 
on group manifolds. To this end, it is crucial for the gauge to be "local" , in the sense that 
the inversion formulas should involve only the fields in a neighbourhood of the origin. The 
formulas of this subsection match this condition. In the following chapters we shall see 
that, giving up the regularity of the fields at the origin and exploiting the residual gauge, 
it is possible to write down some formulas which are similar to those of this chapter, but 
contain integrals on the domain (1 — oo). In that case, the "locality" of the gauge condition 
is lost. 

We remind that the Riemann tensor and the torsion tensor are given by the expressions 

R-bnv = d^lv ~ 9 v r^ + r^rgy - r^rg M , (3.15) 

S% = <V" - d v e% + elTl - e\Y a hv . (3.16) 
Multiplying (3.15) by £ M we have from (3.7a) 

e^(o = e M w(o+a(o- 

If we now put £ — > A£, we obtain 

^(Ar^(A0) = A^ 6 V(Ae) 

and integrating we finally have 

n,(0=e [ i^(A£)AdA. (3.17) 
Jo 

In a similar way, by multiplying (3.15) by £ M and taking into account condition (3.7b) we 
obtain 
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by the same procedure we obtain 

^(A(e»(A£) - O) = \eri(XO + XPS^iXZ) 

and 

em = k+ f\enA^)+^s^(xo]xdx. (s.is) 

Jo 

By substituting (3.17) into (3.18) we have 

e a m = 6 a „+ea b [ R a b ^m-^dx+e f s^x^xdx. (3.19) 

Jo ' Jo 

The equations (3.17) and (3.19) are the analogue of (3.2). The gauge potentials they give 
satisfy the gauge conditions (3.7a) and (3.7b) as a consequence of the antisymmetry of 
and S® u with respect to the indices p, v. However, they solve (3.15) and (3.16) only 
if the functions R^^iO an d S^(^) satisfy some conditions. These conditions have been 
derived and used in the Yang-Mills case [Durand e Mendel, 1982] and a similar treatment 
can be given in the case under consideration. If we substitute (3.17) and (3.19) into (3.15), 
(3.16), after a long calculation we find that they are equivalent to the following projected 
Bianchi identities 

£ M (d^Rbvp + ^c^Rbvp ~ ^b^Rcup) = 0' 

{p,up} 

£ M ^2 (^^p + ^bpStp ~ e^Rbvp) = 0, 

{p,up} 

in which the potentials and are replaced by the integral expressions (3.17), (3.19). 
We have indicated by Y2{nv P } the sum over the cyclic permutations of the indices //, u, p. 

These conditions have a non local character, since they are expressed by integro- 
differential equations. Also the Einstein field equations e% R bjlu = contain the potential 
and therefore take a non local character if we want to express them in terms of the 
curvature alone. In conclusion, if we try to use the inversion formulas to eliminate the gauge 
potentials from the theory (like in Mandelstam's formalism), we get very complicated non 
local field equations. 

3.3 Translation of the origin. 

If we adopt the radial gauge, the geometry of the space-time manifold in a neigh- 
bourhood of the origin Po is completely described by the functions e"(£), r^(£), which 
have to satisfy the gauge conditions (3.7a), (3.7b) and the field equations of the theory 
we are considering. In a similar way, the matter fields are completely described by their 
components with respect to the tetrads, expressed as functions of the normal coordinates 
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£. The only arbitrary elements in this description are the choice of the origin Pq and of 
the tetrad e a (Po). Thus the fields in the radial gauge can be considered as "observable" 
as the fields in Minkowski space-time described by their components with respect to a 
given inertial coordinate frame. In fact, also in this case the values of the field components 
depend on the choice of the origin and of the directions of the coordinate axes. 

The transformation properties of the fields with respect to the Poincare group permit 
to compute the components in the new reference frame in terms of the old components. 
For instance, if V a is a vector field, we have 

V' a {x') = [L- 1 ] a b V h (x), x = Lx' + x . (3.20) 

In the following we generalize this formula to a curved space-time with radial gauge, namely 
we derive the explicit form of the transformations of the fields caused by a translation of 
the origin. We consider a radial gauge with origin at the point Po , coordinates £^ and 
tetrads e a (£) and we start from the point Pi with coordinates £f and from the tetrad 

e' a (0) = L b a e b (Ci) 

to build a new radial gauge with coordinates £ M ' and tetrads e' a (£'). If, for simplicity, we 
take L = 1, the coordinates and the tetrads are connected by 

e" = s"(e 1 ,o, e' a (£') = tta(£i, £>&(£)• 

It is easy to see that for general values of L these relations are modified as follows 
The transformation property of a vector field is 

v a \o = [n-^i.LOi- 1 ]?^) (3-21) 

and tensors of arbitrary order transform in a similar way. In Poincare or Euclidean gauge 
theories it is easy to write also the transformation properties of spinor fields. In a flat 
space-time we have 

3"(£i,0 = # + n*(£i,0 = # (3-22) 

and the transformation property (3.21) takes the form (3.20). It follows from the definitions 
that (3.22) holds for general spaces when £1 and £' are proportional. 

The gauge potentials, which describe the geometry, transform in the following way 
(for L = 1): 

e ?(0 = ^Qp^'h n-HtuZ'Wete), (3.23) 

c At (n = [^- 1 (6,n]c 



(0 + 



(3.24) 
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Note that these transformations do not form a group because the quantities £') and 

^a(Ci)C') can depend also on the initial point Pq- They can be computed by means of a 
method similar to that used in Section 3.1. The result are the following integral equations 



Jo 



- t' d t ^(£1, AO<(S(ei, A£'))lt A£'))^(£i, A£') dA. 



Jo 

It is useful to consider an infinitesimal displacement £1 of the origin and put 

s" (6,0 = er + r + ewco + o(^ 2 ), 
^ a (^a = s b a +^B b a ^')+o(e 1 ). 

From eq. (3.22), which holds when £i and £' are proportional, we get the conditions 

e^(o = o, e"sj M (o = o. 

Formulas (3.21), (3.23) and (3.24) take the form 

r£(0-r? M (0 = efW M + o(e?), 

where we have introduced suitable definitions for 5 p V a , 5 p e a ^ S p T^ [M. and Toller, 1990]; 
for instance, 

&pV a = -Bi p V b + W + A v p )^. 

We may easily obtain integral equations for the quantities A and P, which can be 
solved perturbatively for small values of — 5% and T^. The formulas obtained in 
this subsection turn also out to be useful in applying the radial gauge to the problem of 
non-singular homogeneous and isotropic random fields. 

3.4 Solutions of (2+l)-gravity. 

In this section we give a brief account of the applications of the radial gauge to (2+1)- 
gravity due to Menotti and Seminara. 

In the last years the amount of field-theoretical work about 2- and 3-dimensional mod- 
els has greatly increased. The motivations for this interest are various. In general, lower 
dimensional theories are more easily solved in a complete and rigorous mathematical fash- 
ion than higher dimensional ones. Some important and general effects, like spontaneous 
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breaking of gauge symmetry, anomalies, solitons, have first been discovered in two di- 
mensions. For another class of phenomena, the generalization to higher dimensions is very 
difficult or impossible; there exists, however, some real physical system to which the theory 
can be applied. This is the case, for instance, of the high-temperature superconductivity 
and the quantum Hall effect. 

In the case of (2+1) -dimensional gravityf the main interest of the theory resides in 
its connection to the topological theories and to the problem of the cosmic strings in 4 
dimensions [Vilenkin, 1985]. Starting from the work of Deser, Jackiw and 't Hooft [1984] 
a number of papers has been devoted to the classical solutions in (2+1) dimensions, and 
some new solution techniques have been developed [Clement, 1985, 1990; Deser and Jackiw, 
1989; Grignani and Lee, 1989]. 

All the known solutions, and a lot of new ones, have been re-obtained by Menotti and 
Seminara [1991 a, 1992] using the radial gauge. In fact, we recall that in radial gauge the 
tetrads and the connection can be written - in the absence of torsion - as integrals of the 
Riemann's tensor R a * v (eq.s (3.17), (3.19)). Moreover, in (2+1) dimensions the Einstein's 
equations state a linear dependence between the tensor R a ^ v and the energy-momentum 
tensor Tg\. If the energy-momentum tensor is given, we can therefore write the metric as 
an integral of T£. Of course, in order the whole procedure to be consistent, T£ cannot be 
arbitrary; it must satisfy a constraint relation which is related to the Bianchi identities. 
We shall derive it here following a slightly different approach from that of Menotti and 
Seminara. Let us consider the Einstein action in the first-order formalism, in any dimension 
iV > 3: 

C — nrmt; t / rl N R ab p Cl C 2 pCiV-3 p d fj,vp 1 p 2 ...p N -3<T 
° — WJllbl. I U J> t pi c p2 ••• e p JV _ 3 Co- fcabcic 2 ...Cjv_3<i fc J 

where R a ^ v is defined in terms of the connection Y a p like in (3.15). Coupling the vierbein 
to an energy-momentum source T% and minimizing S we obtain Einstein's equations in 
the form 

^jxv^px ■■■ t -abc 1 c 2 ...c N - 3 d t — lAJll&l. ± d , 

or 

Rfiis e P \ ■■■ £abc 1 c 2 ...c N _ 3 d ^ PlP2 PN 3 Cere = COnst. Trf e . (3.25a) 

Now, it is generally believed f that in gravity it is not possible to "give the conserved source 
and solve for the fields and Y a p ". In fact, the conservation of T°[ cannot be checked 
without any knowledge of the geometry. So in general one must consider the simultaneous 

f For an introduction to the classical and first-quantization aspects see Jackiw [1989]; 
for the second-quantization aspects see Witten [1988], and references. 

| Because of this, there is no field propagation in empty space. This shows that (2+1)- 
gravity is physically very different from the "true" (3+l)-dimensional gravity. 

f See, however, the paper of Boulware e Deser [1976]. 
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evolution of the fields and the sources must (ADM formalism). Let us then look at the 
"inverse problem": given the field, we want to find the conserved sources that produce 
it. By virtue of the Bianchi identities, the tensor T% defined by (3.25a) is covariantly 
conserved. We just need to impose the symmetry, obtaining the constraint 

{j^pv^pX ■■■ e abc 1 c 2 ...c N _ 3 d S pU PlP2 ' " pN ~ 3<J e ae ] ^ = 0. (3.256) 

The task of finding solutions of (3.25b) in any dimension is very difficult. Let us take, 
however, N = 3; eq. (3.25b) becomes 

{R*U abc e^e pd } [cd] =0. (3.26) 

We now fix the radial gauge. The most general structure of a radial connection is 

rf(0=e abc e^eA:(0, (3-27) 
where A a c is an arbitrary tensor field. From (3.27) we have 

R a X e ahc = 2A P + £^ A p c - ? (d p A p - \ e abc e^ v C A pa A"' 

We recall (eq. (3.18)) that in radial gauge the vierbein can be expressed as 

e pd (0 = S pd + I eUA£ c rc d (A£). 
Jo 

Then the symmetry condition (3.26) becomes 



-abp 



2 A bp + ^ A bp - ¥ (d»A b » - l - e mnb e aliV T A p m A^j 

+ e pav r (2 A bp + ^ A bp ) f 1 dX A 2 (£ a ^(A£) - £ b A v a ()4)) = 0. (3.28) 

Jo 



This is the constraint equation given by Menotti and Seminara. 

Some solutions of the constraint in (2+1) dimensions. 

Let us first define the scalar operators 

d d 



The operator D n is invertible, for n > 0, provided it acts on functions which are regular 
at the origin. In fact, if D n f = 0, then / is an homogeneous function of degree —n. The 
inverse of D n has the form 

ATV(0= / e/A A n-1 /(A£). (3.29) 
Jo 
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We shall consider two simple "Ansaetze" for A a ^ [Menotti and Seminara, 1991 a], which 
make vanish the term with the integral of the constraint (3.28) f : 

(1) A%{£) = e wi th UO arbitrary; 

(2) A^(0 = x M/ 9 a (0, with p a (0 arbitrary. 

Substituting the first ansatz into (3.28) we obtain 

D±r(0=d P F(0, (3.30) 

where F is any function. From (3.29), (3.30) it follows that / p (£) is proportional to £ p ; we 
then conclude from (3.27) that the connection is identically zero. So the first Ansatz has 
to be rejected. 

Substituting the second Ansatz, we get the equation 

X " D 2 p a - D x p a = X a D 2 f - e D x pP. (3.31) 

Multiplying both sides of (3.31) by Xn we obtain (if x 2 7^ 0) the vector equation 

X 2 D 2 - (xl) D x ] p = xD 2 (xp) - £D x (xp). 

In a suitable coordinate system, the operator in square brackets becomes an operator of 
the form D 2 , and is thus invertible; it is sufficient to choose as new coordinates z\ = (x£)> 
Z2 = z 3 = (^0? being a and b two vectors which are orthogonal to x an d between 

themselves. We then obtain for p the structure 

P = fx + gl (3-32) 

It can be shown [Menotti and Seminara, 1991 a] that p has the same structure also in the 
case x 2 = 0. Substituting (3.32) in (3.31) we obtain the following differential equation for 
the functions / and g: 

D A g = -D x f (3.33) 
A simple and interesting solution of (3.33) can be found by setting 

9 = 0; f = f[(S0,(H) 

where a and b are two vectors orthogonal to x- It can be easily shown that the correspond- 
ing vierbein is given by 

<(o = s; + \t 2 e ahc e atlv 6 e xc x d 



f In the absence of this term, the holonomic index of T is lowered, in practice, just with 
5® instead of the whole vierbein; this could be viewed as a linear approximation. In fact, 
the other quadratic term in A vanishes too. 
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where 

0(0 = / 1 rfAA(l-A)/(A0- 
Jo 

From this, through (3.3), we obtain the metric 

g^io = + p <w e aPv X p r x a ho x 

{l + \f (x 7 Xt) - (C 7 X 7 ) 2 ] 0(0 } • (3-34) 

It is easy to verify that x is a Killing vector for the metric (3.34). We thus go over 
to the reference system where x assumes its simplest form. Let us consider the cases of 
time-like x an d space-like x- (For the case of light-like x see [Menotti and Seminara, 1991 
a].) 

(1) x time-like. Setting x = (1, 0, 0) the metric becomes 

#00 = 1; 9oi = 0; 



9ij (0 = ± e & HO (1-4 + G) HO ) , 



where the + sign holds for the diagonal and the - sign for the off-diagonal components. 
Going over to polar coordinates we obtain 



= diag < 1, -1, 



1 - ]- l 2 r 2 Hr, 9) 



1 2' 



where 



Hr,e) 



d\X(l- A) f(Xr cos 0, Arsinfl). 



(3.35) 



Eq. (3.35) generalizes the usual conical metric to an extended distribution of matter 
which does not possess axial symmetry. Defining the defect angle by the ratio between the 
radius and the circumference, we obtain 

A6> = -£ 2 M, 
2 

where 



/ d9 I dr r f(r cos 6>, r sin6>) = / d 2 x ^fgT Q Q. 
Jo Jo J 



M 
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Suitably choosing / we obtain the static point-like solutions of Deser, Jackiw and 't 
Hooft [1984] and the string solutions of Grignani and Lee [1989], Clement [1985, 1990], 
Deser and Jackiw [1989]. 

(2) x space-like. Setting x = (0, 0, 1) and performing the hyperbolic transformation 

^ScoshT; £ = SsinhT, 



the metric takes the form 



g^ = diag < =f5 s 



lT^t 2 S 2 $(S,T) 



-1, -1 



where the - sign holds inside the light cone, and the + sign outside. This is a generalization 
of Rindler's metric. 

The formulas we have written in this section are just an example of the solutions 
that can be found using the radial gauge. Moreover, it is possible to define a "reduced 
radial gauge" [Menotti and Seminara, 1991 a], which is useful in the analysis of stationary 
problems and allows a complete solution of the problem of timelike closed curves in (2+1)- 
gravity [Menotti and Seminara, 1993]. It is even possible to solve the constraint equation 
(3.28) in general form [Menotti and Seminara, 1992]. 



3.5 Motion of test particles in a fluctuating field. 



In this section we shall present another application of the radial gauge, namely to the 
case of test particles moving in a weak quantized gravitational field (in four dimensions). 

We start recalling that, from an operational point of view, the definition of Minkowski 
space-time is based on the possibility of building up an orthogonal network of rods and 
clocks. If a gravitational field is present this possibility is lost, but the equivalence between 
inertial and gravitational mass still allows a geometrical formulation of the theory. The 
idea that fields and geometry are intimately related has proved to be one of the most useful 
and fruitful concepts of physics. We shall thus assume that in the presence of a classical 
gravitational field space-time is described by a differentiable manifold M, endowed with a 
metric g. 

We recall that it is possible to define on M a "geodesic structure" [see for ex. Kobayashi 
and Nomizu, 1969]. This structure is in some sense the physical manifestation of the 
geometrical framework of the theory: namely it is invariant with respect to coordinate 
transformations, and it can be operationally tested by observing the trajectories of many 
free-falling test particles. Owing to the equivalence principle, the mass of such particles has 
no importance, as long as they do not perturbate the field. The motion of test particles in 
an external field has been studied in its various aspects by many authors [see for instance 
Papapetrou, 1951; Souriau, 1974; Toller, 1983; Toller and Vaia, 1984]. 

In this section we want to show - as an application of the radial gauge - that a 
weak quantized gravitational field introduces a difficulty in the definition of the geodesic 
structure of space-time at very small distances: namely the vacuum correlations of the 
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field will influence the motion of a test particle, and this influence will depend on the size 
of the particle itself. Such a phenomenon is not difficult to be intuitively understood. In 
fact, the size L of a particle is most properly defined by diffraction experiments. This 
means that when the particle travels in a fluctuating gravitational field its motion will be 
affected just by the fluctuations of wavelenght A such that A > L. The point we would 
like to clarify is the average effect of these fluctuations on the geodesic structure measured 
by the particle. The conclusion is that such an effect is non vanishing and proportional 
to £ pianck /L [M., 1992c]; this confirms that a geometric theory based on the equivalence 
principle is not operationally meaningful at distances comparable with £pi anc k, since the 
geometry depends on the size of the test particles. 

The calculation relies on eq. (3.19), which gives the tetrad e^(£) in radial gauge as an 
integral of R^^iO ( we assume that no torsion is present). Writing -R^(0 m t erms of the 
Riemann's tensor through (3.6), we obtain the equation 

e a M =51 + ^ b f ds s(l - s)R^M) 
Jo 

This is an integral equation for e"(£), which can be iteratively solved to the desired order 
in R. After two iterations we find 

= ^ + K Ft* f ds 8 {i - s)R^(st)+ 

Jo 

+^re /3 er/ ^ 3 (i-*) / ^(i-^^^o^^o + o^ 3 ) 

Jo Jo 

(3.36) 

The crucial observation now is that the Riemann tensor of a weak euclidean field 
9nu = 5 nv + h uv is invariant with respect to gauge transformations. Thus we can substitute 
in (3.36) the Riemann tensor R(£) in the radial gauge with that in the Feynman-De Witt 
gauge, say R F (x). To this end we also must express £ in terms of x, but this will just affect 
the term with one single R. The function x M (£) is given by (3.13). Finally, eq. (3.36) can 
be rewritten as 

<(o = ^ + k ee f ds 8(i - s)R F : 0fl (so+ 

Jo 

-[ dss 3 (l-s) [ dt(l-t)fZ(s,t,0 + o(R 3 ), (3.37) 
jo Jo 



where 



t, = 5i r^er R F 1^(x) r F lM - R F l^ x ) <l» , (3.38) 

L J x=st,, y=st^ 



Eq. (3.37) describes the geodesic motion of the parallel-transported tetrad. Note that 
e^(£) still satisfies the gauge condition (3.7b). Let us now suppose that the gravitational 
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field consists of weak fluctuations quantized around a flat background. We can find the 
"average motion" of the tetrad by averaging (3.37) on the vacuum state. To this end 
we just need to replace the quantities in square brackets with their Feynman-De Witt 
propagators. The term linear in R will vanish in this approximation. One thus starts from 
the propagator of the metric in dimension 4 [De Witt, 1967 c; Veltman, 1976] 

where k = VlOnG is the Planck length and 

P^iupa = ^iip&vo ^iio&vp $p,v$po- (3.40) 

Then one exploits the linearized expressions for the connection and the curvature 

= \& af *{d v Kp + d^Kp - dphpv), (3.41) 
Refill = ^^(d^d^hpv - dpd^h lv - d^d v h^ + dpd v h 1[x ) (3.42) 



and the formula 



esfdadfadv^ = ^(4^ - e <w (3.43) 

In this way one finds, after a long but straightforward calculation 

(fVstO) =— - (344) 

As we explained above, we must now take care of the finite dimension of the tetrad. It 
is apparent that (3.44) will give a divergent contribution after integration in s and t. This 
happens because the arguments of the fields in (3.38) coincide in some points. However, 
if we cut off the modes of the field with wavelenght smaller than L, the propagator (3.39) 
can be replaced by 

(h F {x)h F {y)) ~ {x _y ? + L ^ (3-45) 

we implement this condition restricting the integration in s and t to a domain such that 

\st-stt\>L, (3.46) 



that is 



s(l-t)>|C|- 1 , where C = | > 1- 
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Executing the integration we finally find 



«(C)) = + ^ 0(|C|) - ^) + o(£ 2 ), (3.47) 

where <^(|C|) is an analytic function which vanishes at C = 1 an d whose form depends on 
the details of the regulator. Note that ( expresses the distance in the new unit length L. 
In order our procedure to be meaningful, ( must be not too large. 

Equation (3.47) is the main result. It shows that the "physical" tetrads, parallel 
transported in the vacuum, differ from those of flat space by a correction proportional 
to the square of £pi an ck/L. Our position has been, of course, to consider the Feynman- 
De Witt gauge just as a way to fix the dynamical components of the field, whereas the 
true coordinates have to be constructed in an operational way. If one accepts this point 
of view, one concludes from eq. (3.47) that at very small distances, when L becomes 
comparable with £, the size of the test particle influences its average motion in a fluctuating 
gravitational field. 
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4. RADIAL PROPAGATORS. 



In this chapter we shall introduce a technique [Menotti and Seminara, 1991 b; Menotti, 
M. and Seminara, 1993] which allows to obtain the propagator of a gauge field in a ar- 
bitrary "sharp" gauge (that is, a gauge fixing obtained by insertion of a delta function 
in the functional integral), starting from the Feynman gauge. This technique, called the 
"projectors method" is very general, since it works in arbitrary dimension and can be 
easily specialized to electrodynamics, Yang-Mills theory and Einstein gravity in the first- 
or second-order formalism. It also allows to select the various solutions of the propagator 
equation, which are usually connected by a residual gauge transformation and have differ- 
ent regularity properties. In the case of the radial gauge, as we shall see, there are three 
propagators, which differ by the behaviour at the origin and at infinity. In Section 4.1 we 
shall explain the method in generic form, with emphasis on the algebraic properties; in the 
Section 4.2 we shall apply it to electrodynamics, specifying also in details for illustrative 
purposes the convergence and regularity properties, and in Section 4.3 we shall apply it 
to gravity. Finally, in Section 4.4, after comparing the radial gauge vacuum correlations 
with those of the Mandelstam covariant, we shall compute explicitly their components in 
a simplified case, in order to distinguish the components which are physically meaningful. 

4.1 The gauge projectors method. 

Let us consider a generic gauge field A(x). A can be an electromagnetic field, or a 
Yang-Mills field, or finally the gravitational field, in the first or second order formalism. 
Let the linearized equation of motion for A(x) in the presence of an external source be 
written in the form 

K x A{x) = -J(x), 

where K is a linear, non-invertible, hermitean "kinetic" operator and J(x) is an external 
source coupled to A. The gauge transformations of A have the form 

A(x) — > A(x) + C x f(x), f(x) any function. 

Here C x is another linear operator, which has the properties 

K x C x = (gauge invariance) 

and 

C x K x = =>- C x J(x) = (source conservation). (4.1) 

We assume that we can always add to the kinetic operator K an operator K T which 
makes K invertible, as it happens, for example, in the Feynman gauge. We assume 
to be of the form T ', as deriving from a quadratic gauge fixing J dxJ r2 {A(x)). T is 
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meant to be a linear operator on A. The propagator G T corresponding to this gauge 
fixing satisfies 

(K x + K£) G^ix, y) = -8(x - y) (4.2) 
and has the following property 

j dy G T (x, y) J(y) = if Cj J(y) = 0. (4.3) 

In other words, vanishes when applied to the fields generated by physical sources. In 
fact, applying to (4.2) we get 

Ct K x G r {x, y) + Ct Tl T x G T {x, y) = C\ 8{x - y). 

But using (4.1) we also have 

Cl^ x G r (x,y) = Cl6(x-y) 
and integrating on a conserved source J we obtain 

J dy Cl Tl T x G T (x, y) J{y) = if Cj J(y) = 0. (4.4) 

We notice that T C is the kinetic ghost operator and as such invertible. Then from (4.4) 
we get 

J dyT x G T {x,y) J(y) = 

and multiplying by T\ we finally prove (4.3). 

Let us now impose on A a generic "sharp" gauge condition Q 

A G (x) = {A(x) : <3{A{x)) = 0}. (4.5) 

Q is meant to be a linear function of A and of its derivatives. The field A G (x) can be 
obtained from a generic field A(x) through a (generally non-local) projector P g 

A g (x) =P G [A(x)] = A(x)+C x F g [A(x)}. (4.6) 

We require this projector to be insensitive to any "previous gauge" of the field, namely to 
satisfy 

P G [C x f(x)] =0, for any /(x). (4.7) 
We shall now prove the following properties of the adjoint projector . 
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uces conserved sources. 
For, integrating (4.7) on a current J(x) we have 

J dxP G [C x f(x)}J(x) = 0; 

by definition of the adjoint projector, this means that 

J dx(C x f(x))P g \j(x)]=0 

and integrating by parts we have 

dxf(x)ClP G \j{x)] = 



or, due to the abitrariness of f(x) 

ClP G \j{x)} = for any J(x). 
(2) P gj leaves a conserved source unchanged. 

Let us suppose that J is conserved. Using (4.6) we have for a generic A 

J dxA(x) P g] [J{x)] = J dxA(x)J(x) + J dx {C x F G [A(x)]} J(x). 

Integrating by parts the second term on the r.h.s. and using (4.1) and the arbitrariness of 
A we have 

P g \j{x)] = J(x) if Cl J(x) = 0. q.e.d. 



The equation of motion obtained varying the action under the constraint (4.5) is 

P g \k x A g {x)} = -P G \j{x)}. 
From (4.1) and Property 2 we have 

K x A G (x) = -P Gt [J(x)], 

or for the propagator 

K x G G (x,y) = -P G \d(x-y)], (4.8) 
where the meaning of the r.h.s. is 

J dyP Gt [5(x-y)] J{y)=P G \j{x)]. 
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Next we show that a solution of (4.8) is 

G G (x,y) = (P G [A^(x)]P G [A^(y)]) Q7 (4.9) 
where A T denotes the field in the original gauge. Integrating on a source J(y) we have 

J dyK x {P G [A^{x)]P G [A^{y)}) Q J(y) = 

= J dyK x {P G [A :F (x)}A^{y)) P G \j{y)} = 

= J dyK x ({A^(x)+C x F G [A :F (x)]}A^(y)) Q P G \j(y)} = 

= J dyK x (A^(x)A^(y)) Q P G \j(y)] = 

= J dy(K x + K^)G :F (x,y)P Gt [J(y)}- J dy K T X G T \x,y) P G \j{y)] = 

= - J dyd(x-y)P G \j(y)]. (4.10) 

In the last step we have used (4.3) and Property 1. Note however that G G defined in (4.9) 
remains unchanged if we replace the original field with any other gauge equivalent field. 

Given two different projectors Pi and P2 which project on the same gauge (which as 
a rule differ for different boundary conditions), one has 

P 1 P 2 = P 1 ; P 2 P 1 =P 2 , (4.11) 

due to (4.7). Thus also P12 = a Pi + (1 — a)P 2 is a projector on the considered gauge and 
one can write down the Pi2-projected Green function equation (we omit the suffix Q) 

K x G(x,y)=Pl 2 [5(x-y)}- (4-12) 
It is immediate to verify that a solution of (4.12) is also given by 

a (P 2 [A(x)] Pi[A(y)]) Q + (1 - a) (Pi[A(x)] P 2 [A(y)]) . (4.13) 
Namely, repeating the same procedure of eq. (4.10) we have 

K x a (P 2 [A(x)]Pi[A(y)]) = aK x (A(x) Pi[A(y)}) = a P}[S(x - y)). 

Acting similarly with the (1 — a) term in (4.13), we get (4.12). We notice that for a = \, 
(4.13) is symmetric in the exchange of the field arguments. 

We close this section writing in explicit form the operators which appear in electro- 
dynamics and linearized Einstein theory. 

Electrodynamics and linearized Yang-Mills theory. 
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This is the most simple case. We have the following identifications 



^4 — > A M ; 
J > Jjj,] 
Cf - dp /; 

where is the operator which is produced by the usual Feynman gauge fixing ^(d^Ap) 2 . 

Linearized Einstein gravity in the second-order formalism. 
In this case we have 

A > hp U ; 
T — > T • 

Cf —> (6 a(7 dp + dapdcr) f a ; 

K -> Kp Vpa = ^ [2 5^ M 0„ + 2 ^ <9p<9 CT + 

- (<W + S vp dpd a + 5p, a d v d p + 5 VU d fl d p ) + 
+ (Sp P + <W b~ up - 2 b~p U 5 pa ) d ]; 

k t -> = 1 [-(2 v + 2 <^ 

+ (5 MP d v d a + 5 vp dpd a + 5p a d v d p + 8 va d^d p )}. 



Here K T is the operator produced by the harmonic gauge fixing 

o i ()''h..,. - -r)"l>i' \ 



I^V-^) • (4-14) 



Linearized Einstein gravity in the first-order formalism. 
The quadratic part of the lagrangian has the form 

l {2) = -(C7 dp rf r 7 c + ^ r°„ i? + t* r- + r« 6 ), 

where 

r a = e a - <5 a 

Ct = <W - 

Co 7 = - « + W 
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and and are the energy-momentum source and the spin-torsion source, respectively. 
The gauge transformations have the form 



Cf 



r) ab 8 

u °cd U V 









I Qcd J 



A -<Wd 

The field equations are given by 

KA=J " . . 




a 6 



and the gauge-fixing term has the form 





K T A 



.0 AK^ hv 5^ + 2(3{5 ah 5l-5» h 8Z) 




and is produced by the harmonic gauge fixing (4.14) with h^ v = T^5 av + T®5 a ^, to which 
the symmetric gauge fixing ^(r^5 a u — T®5 a ^) 2 has been added. 

4.2 Electrodynamics. Regularity properties. 

In this Section we shall apply the projectors method to the case of electrodynamics 
or linearized Yang-Mills theory in radial gauge [Menotti and Seminara, 1991 b]. We start 
from the inversion formula (3.2), writing it in the form 

A° l (x) = P°[A} fl (x) = x p [ dX X F p ^(Xx). (4.15) 

Jo 

Let be r = The integral in (4.15) converges if |-F(;r)| < r~ 2+s as r — > 0. We notice 
that if the field A^(x) is such that < r~ 1+e for r — > 0, then (4.15) can be rewritten 

as 

Al{x) = P°[AUx) = A^x) J\\ Xp Ap(Xx). (4.16) 

This shows that the general form (4.6) of the propagator is maintained, provided the fields 
are not too much singular. There is, however, a limiting case in which (4.16) is not true 
any more: namely, when we start from a field which is already radial and differs from 
A {x) by a residual gauge transformation of the form (3.1). For instance, let us take as 
starting field the radial field A^(x) = —x p f x dX XF pp (Xx). Suppose that F pu (x) is finite 
at the origin. It is easy to see that | | ~ r -1 as r — > 0, while F keeps finite. This 
happens because the singular part of A is a pure gauge, which does not contributes to 
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F. The exact significance of this detail will become more clear as we proceed with our 
discussion. In practice, we assume the following definition of P°. 

Definition of P° : 

Al(x) = P°[A]p(x) = x p f 1 dXXF p ^(Xx). (4.17) 

Jo 

In a similar way, we can define a projector in which the integration limits and oo are 
exchanged. 

Definition of P°° : 

/oo 
dXXF P n(Xx). (4.18) 

It is easy to verify that the difference between A and A°° is a residual gauge term of 
the form (3.1). 

It is apparent from (4.17), (4.18) that P° and P°° are projection operators, namely 
we have (P ) 2 = P°, (_p°°) 2 = p°°. In addition, the following properties hold (compare 
(4.11)) 

pO poo _ pO. poo pO _ poo 

The definition domains of P° and P°° are respectively given by 
{F : \F(x)\ < r~ 2+s for r — > 0} and {F : \F(x)\ < r~ 2 ~ s for r — > oo}. Restricting such 
domains respectively to 

{F : \F(x)\ < r~ 1+e for r — > 0} and {F : \F(x)\ < r~ 3 ~ e for r — > oo}, 

we find that the fields A and A°° can be characterized in a simple way through their 

asymptotic behaviour. Namely, we find that 

- A^(x) vanishes at the origin and decreases like r _1 at infinity; 

- A^(x) grows like r~ l at the origin and vanishes like r~ 1_e at infinity. 

In fact, let x be the unit versor of x. We have 

Al(x)=x p [ dX X F PfJi (Xrx) = -f / dttF PfI (tx). (4.19) 
Jo r Jo 

In the limit r-^0we obtain (in the domain above) 

A°(x) <^[t 1+£ ] r ^0 forr^O. 
On the other hand, in the limit r — > oo, the integral converges and A has the form 

Al{x) = - r H%x) forr^oo, 
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where H®(x) is a homogeneous function of degree 0. Analogous formulas hold for A°° . In 
the limit r — > 0, A 00 is regular, in the sense that it vanishes faster than r _1 . At the origin 
we have instead 

A™(x) = --H°(x) forr^O. (4.20) 

From eq.s (4.19) - (4.20) we see that if we compute a radial Wilson loop (fig. 1) going 
to infinity, we obtain the same result both using A and A°°. In the first case the effective 
contribution is localized at infinity, in the second case at the origin. Finally, if we define 
the mixed projector 

P s = I(p° + J p°°), 

we find that the contribution of the field A s to the radial loop is splitted into two parts: 
one half at the origin, one half at infinity. 

For the computation of the adjoint projectors P°\ p°°^ anc [ p s ^ we refer to the 
mentioned work of Menotti and Seminara. 

By means of the projectors P°, P s we can now write down three symmetric radial 
propagators, denoted by G°, G°° and G s : 



G? a ,{x,y) = (P D [A] ll {x)P [A] v {y)) = 

= X Py° [ dXX [ dTT(F p ^Xx)F au (ry)) ; 
Jo Jo 



G™(x,y) = (P°°[AUx)P°°[A] l/ (y)) = 

dXX J dTr(F ptl (Xx)F aiy (ry)) ; 



G%{x, y) = \ (P°[AUx) P°°[A]M) + \ (P°°[AMx) P°[A] v (y)) = 
= -\x p y a J dXxj^ drT(F pp (Xx)F ai/ (ry)) 
-^Y| dXXJ dTT(F pp (Xx)F ai/ (ry)) . 

To compute them explicitly we make use of the formula 

x p y° (F ptl {x) F au (y)) = x p y a D ppua V(x - y), 

where 

J-Juovo — O uo ^ , . n _ ~T Ova ^ .. n „ O ua O, 



' '"' T ilp dx v dy° " n i).ri' On' 1 '"" <).,■'■ i)^ "''().>■" Oif 
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and V(x — y) is the Feynman propagator. It can be shown that the integrals of G° converge 
in dimension N < 4, while those of G°° converge for N > 4; finally, the integrals of G s 
converge for N > 3. Let us compute explicitly, for instance, G s . Using the formula 



/ dp jf — = f 

Jo ' [(px - y) 2 l~ _1 Jo 



dp 



1 



[(px - y) 2 } 2 1 Jo [(px - y) 2 ] 2 
we can write the integral of interest in the form 

r(i-i) i ' < ! 



JV-4 



dp 2Z~ 

o {(x-py) 2 }2 



D(x,y) 



An N / 2 N-A | ,/ 
1 



[(x - py) 2 ] 2 1 



dpp 



dp 



N-4 



[(px - y) 2 ] 2 1 



o ' [(px - y) 2 } 2 1 



(4.21) 



The first integral which appears in (4.21) can be rewritten as follows 



h 



,iV-4 



o dp [p 2 X 2 -2pXYcos6 + Y 2 ] N / 2 - 1 ' 



where X = \x\, y = y , and performing the change of variables p = y ^ , with <i 

ct 

c 1 



VX 2 + Y 2 - 2XY cos 9 , we obtain [Erdely et al., 1953] 

-N-4 



Y^d 3 -" / d( 
Jo 



[e + 2c ( 



) + i] 



N/2-1 



N-3 

being cos (j) x = Y ~ X cose 



1 3-iV / /V — 1 \ 3-JV 

^ y-^ 3 "^ (4sin 2 ^)^ r ^— - (cos^) , 



Here 6* is the angle between x and y, and <p x is the angle opposite 

3-JV 



to x in the euclidean triangle of sides x, y, d. Expressing P 3 2 N in terms of hypergeometric 

2 

functions, we have 

1 d 3 ~ N 



3-N 
2 



2 Vx\ ~ „ , N -3 5-N N -1 . 2 (f> x 

■ cos — 2 F\ , ; ; sin — 

N -3 y- 1 V 2/ V 2 ' 2 ' 2 ' 2 



In a similar way we can compute 



h = 



dp 



o r [p 2 X 2 -2pXYcose + Y 2 ] N / 2 - 1 ' 



obtaining 



— N 



3-JV 

JV-3 X-i V Sm i) 2 ^l^'^^ 003 ^ 



y 3 



iV-3 5-iV AT-1 
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The behaviour of [D(x, y) +D(y, x)] as X — > 0, and y and 6> are constant, is obtained from 
I\ and I2 by noticing that, in this limit, 4> x — > and y — > 7r — 6*. We then have 

3-JV 

1 f 2 #\~ „ 5-iV iV-1 2 0\ 

"iV^( Sm 2 J ^{——;—^ OS 2) X 

This is also the behaviour for x and 9 fixed, while Y — > 00. Owing to the symmetry of 
[■D(x, j/) + -D(y, x)], we have that for Y — > 0, x and 6* fixed 

1 • 2#\^ 5-iV AT-1 2 0\ 

~ sin - 2 Fi , ; ; cos - x 

N-3 \ 2J V 2 2 2 2 / 

( Y 3 ~ N X 3 ~ N \ 

yx^ 1 y^J 

and the same holds for X — > 00, y and # fixed. Hence we see that the behaviour for 
X — > 00 and X — > is the same, except for the sign. Being D^ vpcr a zero degree operator, 
this also holds for the whole propagator G^ v (x, y). 

4.3 Gravitational propagators. 

In Section 4.1 we have developed in general form the projectors method, which pro- 
vides an algebraic way for constructing the propagator in any "sharp" gauge, provided it 
is known in a Feynman-like gauge. Then in Section 4.2 the projectors and propagators 
of electrodynamics were explicitly written, taking care of their existence and regularity 
properties. 

The reason for which it is desirable to impose the radial gauge in gravity as a sharp 
gauge (that is, to insert a delta into the functional integral) is that the geometrical meaning 
of the radial gauge can be formally preserved, at the quantum level, only if the fields 
admitted in the functional integration satisfy the gauge condition. 

The functional integral of gravity in the first order formalism is given byf 

The gauge fixing term is 8{x^T a J } )8{x^{e a p — 8%)). Note that it preserves the symmetry 
between vierbein and connection. 

f It is well known that the rigorous definition of these integrals is affected by the un- 
boundedness of the euclidean action and by the arbitrariness of the functional integration 
measure [see Menotti, 1990]. We shall not dwell on such problems here. 
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As it happens in Yang-Mills theory, the ghosts associated to the local Lorentz symme- 
try formally decouple, while those associated to the diffeomorphisms survive (in the first- 
as well as in the second order formalism). 

We remind that in the functional integral approach the correlation functions are com- 
puted by averaging the products of fields like 0(x)0(y) on all geometries weighted by the 
exponential of the gravitational action. Fixing the gauge to the radial gauge gives x, y,... 
a well defined meaning, as the points that acquire geodesic coordinates x, y,... in each of 
the geometries we are summing over. Hence, when the points lie on a straight line through 
the origin, their correlation is automatically geodesic. 

Next we write the expressions for the radial projectors which are the analogue of the 
P°, P°° and P s of the preceding Section. 

P° is given by (3.17) and (3.19) substituting to R a ^ u and their linearized expres- 
sions 



R 

S 



L a 



(x) 



d u T ah 



= d^(x) - d u r;(x) + r a »{x)8 bv - Vl\x)8 h ^ 
which are invariant under linearized gauge transformations. We thus have 

/ 



V r«(x) 7 



x v \ dXXR^Xx) 
J o 

v x b I dXX{l-X)R^{Xx) + x" [ dXXS^(Xx) 
Jo ' Jo 



\ 



(4.22) 



V r«(x) / 



x" I dXXR^(Xx) 



/oo roo 
dX A(l - A) R^(Xx) - x u J dX X S^(Xx) ^ 



(4.23) 



The projector P° produces radial fields which are regular at the origin (i.e. behaving like 
the original fields) and gives a connection r^ 6 (x) behaving like 1/r at infinity. On the 
other hand P°° projects on a radial field which is regular at infinity and gives a connection 
behaving like 1/r at the origin. 

P s , which allows us to construct a finite propagator in iV > 2, treats the origin and 
infinity in symmetrical way by giving the same 1/r behavior in the two limits with opposite 
coefficients. This can be shown in a similar way as we did for electrodynamics [Menotti, 
M. and Seminara, 1993]. 

The adjoint projectors can also be easily computed and their expressions can be found 
in [Menotti, M. and Seminara, 1993], like many other details that we shall omit in the 
following. 

When we construct by these projectors the propagators 

(P°(r,r)P°(r,r))o 
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and 



(P oo (r,r)P oo (r,r)) 

we find that the first one is always divergent, while the second one diverges for N < 4. Thus 
we shall construct the solution for the P s projected Green's function equation, namely 

8&Z*d x ) V G d ^ 



t / W , 

= -P s1 \ \S N (x-y). (4.24) 
by use of the now familiar procedure. One shows that the propagator 

i/po p°°(rg;„ ,(y),T«;(y))) + (P o ^p° o ) (4.25) 



2 



r M a (x) 



is a convergent symmetric radial solution of (4.24) for all N > 2. The explicit form 
of solution (4.25) of the Green' s function equation (4.24) is easily computed by using 
(4.22) and (4.23) where the correlators between Riemann and torsion two-forms, which 
are invariant under linearized gauge transformations, can be obtained using e.g. the usual 
symmetric harmonic gauge. 

Let be f (x, y) the ultra-local part of the propagator of the anholonomic connec- 
tion (that is, the part which is proportional to 5 N (x — y)). M is given by the expression 
[Menotti and Pelissetto, 1987] 

M a ^ cd (x, y) = - % - (V 5^ d + 8<* 5% - ^_ 8* Sf^j 5 N (x - y). 

(8 ah ' cd is an antisymmetric symbol in the pairs (a, 6), (c,d).) 
We have 

(R L t(x) R LC P d M) = (R L t(*) R L t{y)) n + [d,d p M^ d ( x , y )} ^ ] [pa] ; (4.26) 



(s L lM s Lb P M) = lKT^y)^^\ uUpaV (4.27) 



(R L t ( x ) SL lM) = %M^ cd {x, y) S d<7 ] [H [pCT] . (4.28) 
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In (4.26), ( ) TI denotes the correlator in the second order formalism: 

(R»vap(x) Rpa\y(y)) = X 



rpV ra'8' cX'8' . c\'u' so? 8' X p'8' 2 raVwV'rAV 

5 pa KB S X-y + S Xj °ad 5 pa ~ ^ S aB Ka 6 Xj 



X 



d^d a >dp>d\>V(x - y). 

The ultra-local nature of the correlators (4.27) and (4.28) reflects the well known fact 
that the torsion does not propagate in the Einstein-Cartan theory. 

The integrals over A and r of the (RR) 11 correlators are similar to those given for 
electrodynamics (Section 4.2). The contact terms generated by M give rise to integrals in 
A and r which are convergent. In fact the generic integral is of the form 

pl poo 

J d\X A J dTT B 5 N (\x-Ty) = 

pi poo 

= J d\X A J dTT B 5(\\x\-T\y\)(\\x\)- N+1 5 N - 1 (n x -n y ) = 

= s N -\n x - n y )e(\x\ - \y\) A + B 1 _ N + 2 * 

{lyl-^lxf-^ 1 - |y| j4 " JV+1 |a;r A " 1 ) • 

We notice that in three dimensions the correlator (4.26) is identically zero, as can be 
explicitly verified by observing that the matrix component (2,1) of eq. (4.24) takes the 
form (R L (x)r(y)) = 0, which implies (R L (x)R L (y)) = 0. This means in turn, through 
(4.22) and (4.23), that (T(x)T(y)) = 0. Hence (r(x)r(y)) vanishes, except for singular 
contributions which arise when the origin, x and y are collinear; such contributions are 
given by (4.26), (4.27) and (4.28), that is, by the torsion- Riemann and torsion-torsion 
correlators. 

In the second order formalism, the Riemann-Riemann correlator, which in Minkowski 
space corresponds to the T* product, is not identically zero even in dimension N = 3. 
Namely, if this were the case, all propagators would vanish identically; but this cannot be 
true, since the collinear singularity must be present in the propagator (h(x)h(y)) in order 
to produce the conical defect which is typical of 3D gravity. 

A merit of the radial gauge is to show explicitly the absence of propagation in the 
three-dimensional theory. The harmonic gauge - on the contrary - propagates a pure 
gauge field. 

4.4 Relation to the Mandelstam covariant. Significant components. 

We recall that the "Mandelstam covariant" TZ a bcd{x, P) (Section 2.1) is the curvature 
tensor observed in the local reference frame parallel transported to a point from a fixed 
origin x along the path P. In the version of Tsamis and Woodard (Section 2.2) the path 
P is assumed to be geodesic. 
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It is then clear that in the radial gauge the (geodesic) Mandelstam covariant referred 
to the origin of the coordinates can be written as 



where R a b P v is given by eq. (3.15). One verifies from this expression that 1Z is a scalar 
under coordinates transformations. 

Thus the (geodesic) Mandelstam covariant is a composite field which can be con- 
structed by radial fields. Its correlations can be written using the propagators of e^(£) 
and r^ 6 (£) given in Section 4.3, and the resulting analytical structure will be the same, 
including hypergeometric functions etc. The index structure is quite complicated how- 
ever, and it is more interesting instead to show explicitly the components of the radial 
vierbein-vierbein correlation function, in order to gain a feeling of its behavior. 

Denoting as usual the vierbein by e a ^ = 5 afl + r a/J and all the integration procedure 
of the propagator (see Section 4.3) by J d[X] J d[t], we may write 



(r a M) Tc P (y)) = / d[X] / d[t] [x»x b y°y d (Ra b »Ax)R cdp Ay))} x=Xx . y=ty • (4-29) 



In this formula, R a b^( x ) denotes the usual "holonomic" components of the (linearized) 
Riemann tensor, that is, R a bnu(x) = $a$b R a p P v(x), and similarly for R c d P a- Also, "ul- 
tralocal" correlations are disregarded (compare eq. (4.26)). 

According to our usual assumptions, the two vierbein fields have to be connected by 
a geodesic of length £; this can be realized for instance by choosing the radial coordinates 
x and y in the following way: 



This singles out the direction 1 and causes r to vanish by radiality when one of its indices is 
equal to one; the other three directions, however, are completely equivalent. The expression 
in square brackets of (4.29) reduces then to 



and there are three possibilities of choosing the indices {a//} and {cp} of the two vierbein 
fields: 

(a) like in (T22T22) or (T23T23) ...; this represents the correlation between the same compo- 
nent of the same vector of the vierbein; 

(b) like in (T22T23) or (T22T32) ...; this is a correlation between different components of the 
same vector or between the same component of different vectors of the vierbein; 

(c) like in (T22T33) or (T22T44) ...; that is, a correlation between different components of 
different vectors, but chosen in such a way that a = p and c = p. 





(4.30) 



£ 4 (Ral P l(x Q ) Rdpl(yo)) 



(4.31) 
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It is straightforward to insert these choices of indices in (4.31) and to evaluate the 
propagator for each choice. Also, the integrals in A and t are trivial in this case and the 
dependence on £ factorizes as £ -2 , as expected on dimensional grounds. The three cases (a), 
(b) and (c) give different results, and it is interesting to look at them keeping in mind that 
each vector of the vierbein is parallel-transported between xq and 2/0 in any configuration 
of the functional integral used to produce the averages. It is found the following: 
The correlation (a), which can also be expressed as (T^Ty) (no sum), does not vanish. It is 
in fact the basic correlation function of the radial gauge, as we shall now explain in short. 
Namely, it will be shown in Section 5.3 that the invariant correlation between two vierbein 
fields includes the matrix U of the parallel transport: 

C Vierbein (-D) = (ejW'^.^M) 

(here D is the distance between x and x'). But in radial gauge this matrix reduces to the 
identity along any ray, thus we have for the choice of coordinates (4.30) 

GVierbein(0 = ^M^fe/o)) + (t^^O^^o)) + {T2i(x )T^(y )) + 

+ (733 (£0)733(2/0)) + (7-34(^0) 7-34(2/0)) + (7-44(^0)7-44(2/0))- (4.31) 

This is true to any order. We notice that the terms of (4.31) which differ by the first 
index of the vierbein are certainly equal, since the vectors of the vierbein are equivalent. 
Moreover, if no external source is present, different components of the same vector are 
equivalent too, in the average. One concludes that all the radial correlations of the form 
(rijTij) are equivalent on the vacuum and equal to g of the invariant correlation GVierbein- 

The correlation functions (b), also expressible as {r^rkj) or {r^Tik) (no sum), vanish to 
lowest order and are in general not significant, because there is no direct connection (apart 
from the orthonormality realtions) between the values of two different components of the 
same vierbein vector or the values of the same components of two different vectors. 

The correlations (c), also expressible like (raTjj) (no sum; i ^ j), do not vanish to lowest 
order. Nevertheless, this property is due to the weak-field approximation and does not 
hold in general. It is easy to verify by geometrical arguments that on a almost-flat space 
the i-th component of the i-th vierbein vector keeps constant in sign, so its correlation 
does not vanish. This ceases to be true on a generic strongly curved space. 
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5. GAUGE INVARIANT CORRELATIONS. 



In the first part of this review (Chapters 2-4) we have shown how some correlation 
functions of the gravitational field at geodesic distance can be computed in a special 
physical gauge. By this technique two-point functions are obtained, which really depend 
on geodesic coordinates x and y (more exactly, they are hypergeo metric functions of x and 
y - see Sections 4.2, 4.3) and have "parallel transported" tensorial components which are 
a generalization of the usual Lorentz tensorial components (compare Section 3.3). 

In this second part (Chapters 5, 6), on the contrary, the most relevant scalar, or 
"gauge invariant" correlations will be introduced; by this we mean functionals of the field, 
whose corresponding classical quantities are independent of the choice of the coordinates. 
The most typical and useful example of this kind of correlations is the Wilson loop of the 
Christoffel or anholonomic connection (Sections 5.1, 5.4). 

The physical interpretation of these quantities turns out to be more clear and interest- 
ing than in the case of the gauge-dependent correlations. On the other hand, the "geodesic 
corrections" are less easy to be taken into account in this case; we shall illustrate a special 
procedure for that in Section 6.3. 

This Chapter is organized as follows. In Section 5.1 the matrix U of the parallel 
transport is defined, both in terms of the Christoffel connection T^ v and of the gauge, 
or "anholonomic" connection T a ^ h . In Section 5.2 we specify the dynamical scheme as 
the traditional perturbative scheme which starts from the Einstein lagrangian, splits the 
gravitational field into a flat "background" and a weak, quantized part. In Section 5.3 the 
invariant two-point functions involving the Riemann curvature are computed to leading 
order. They vanish as a consequence of the equations of motion. On the contrary, the 
correlations of the metric and of the vierbein do not vanish. Finally, in Section 5.4 we 
compute to leading order the Wilson loop of the Christoffel connection and find that it 
vanishes too. The physical interpretation of this property is postponed to Chapter 6. 



G.l Geometrical definitions. 

We consider a classical spacetime M described by a metric tensor g^ v {x) of signature 
(— 1, 1, 1, 1) (the conventions are those of Weinberg [1972]). 

The variation of a vector V a by an infinitesimal parallel transport is defined by 

dV a = -T?Jx) V 13 dx^, (5.1) 



where F®p is the Christoffel connection 



Integrating (5.1) we find that the parallel transport of V along a finite differentiable curve 
connecting the points x and x' is performed by the matrix 

^(x,x / )=Pexp I dy»r«p(y), (5.3) 

J X 
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where the symbol P means that the matrices (r M )Jg = are ordered along the path. 
The indices of Up (x, x') are lowered and raised by g ai {x) and g l3 ' y (x'), respectively. 

When the manifold is curved, the matrix U depends not only on the end points x 
and x', but also on the path. So, if C is a smooth closed curve on M, we define the loop 
functional (or "holonomy") W(C) as 

W(C) = -4 + TrU(C) = -4 + TrPexp^ dxT^x). (5.4) 

The term —4 sets the holonomy to zero in the case of a flat space, when the matrix 
U reduces to an identity matrix. 

Under a coordinates transformation x — * £, the matrix U transforms in the following 

way 

U%{x,x') ^U$(x,x') 

For a closed curve, this transformation, being of the form 

U -> OWO" 1 , 

does not affect the trace of U. So the loop W(C) is invariant with respect to coordinate 
transformations. 

In the first order formalism, the "anholonomic" components of a vector are defined 

by 

ya = yv e a( x y 

The equivalent of (5.1) in terms of the anholonomic connection Y a ^ h is 

d ya = _Y^ X )V h dx^ (5.5) 

and the matrix U of the finite parallel transport has an expression which is formally the 
analogue of (5.3), namely 

U^x,x') = Pexp ! dy»r; b (y). 

J X 

We remind (see eq. 3.5) that the relation between the connections and Y a ^ h is the 
following 

L H/3- e a e M f (3b + e a°l3 e v 

and that the relation between the matrices Up and U£ is 

US(x,x') = e a a (x)U^(x,x')e^x'). (5.6) 
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dx c 



dx p 

w 



It is known that gravity in the vierbein formalism has a local Lorentz invariance, since 
the definition of the vierbein, eq. 3.3, is insensitive to a Lorentz rotation of e a (x), e b (x). 
The connection Y a ^ b is then completely analogous to an usual gauge connection, and its 
Wilson loop 

W(C) = -4 + Tr(W 6 a )(C) 

is a natural invariant quantity of the theory. But from eq. (5.6) we see that this loop is 
equal to that defined in (5.4). So the Christoffel connection T"^ and the anholonomic T a ^ h 
connection have the same loop, denoted by W(C). In the computations we shall employ 
the connection T"^, which is usually simpler to deal with. 

When the exponential in (5.4) is expanded, one obtains terms with 1, 2, 3, ... fields 
T. We introduce the notation, to be used in the following 



U = 1 + £ dx» + 1 P j dx» j dy v r M (x) V v {y) 



+ 



and 



= 1+WW + J^ + ... (5.7) 



W= -4 + TrW = TrW (1) + ^TrW (2) + ... (5.8) 



5.2 Dynamics and perturbation scheme. 

Let us now introduce dynamics through the Einstein action, which has the form (for 
any N > 2) 

S =^2 J d N x^gJx)R{x). 

A completely consistent quantum theory of this model does not exist yet. In view of 
applications to lattice gravity, we shall refer in the following, as a possible approximation to 
the full theory, to the (regularized) functional integral approach (see for instance [Hawking, 
1979; Mazur and Mottola, 1990]; compare also Chapter 7). 

In the perturbative evaluations of this chapter, however, we follow the "traditional" 
approach and regard quantum gravity as an ordinary field theory on a fixed flat back- 
ground. The fundamental field is Kh^ u (x) = g^ v {x) — <5 Mi ,, which is subject to the gauge 
transformations 

hpvix) -> hp V {x) + dpUix) + duf^x), (5.9) 

where / M (x) is an arbitrary function. Eq. (5.9) represents the action on h of a linearized 
diffeomorphism x" 1 = / M (x). It should be noticed, nevertheless, that the argument x re- 
mains unchanged in the transformation (5.9): this is what is meant by "fixed background". 
The transformation of V which corresponds to (5.9) is 

r" (*)->r« (x) + d lt d v f a (x). (5.io) 
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Keeping the quadratic part of S and adding to it the harmonic gauge-fixing 

we obtain the Feynman-De Witt propagator (compare Section 3.5, where however h differs 
by a factor k) 

(V(*)VG/)> = -j^y 2 5 ^+^p- 5 ^ . (5.11) 

Due to the non-polynomial character of the lagrangian, there are infinitely many 
interaction vertices; the first two ones, respectively proportional to n and k 2 , connect 3 
and 4 fields h. Hence the first few orders of perturbation theory are formally very similar 
to those of Yang-Mills theory [see for instance Veltman, 1976]. 

As it is known, the Einstein action is not the only action which describes correctly the 
macroscopic behaviour of gravity. In particular, the (R + _R 2 )-action (see Chapter 7) has 
been proposed a long time ago as a renormalizable generalization of General Relativity. 
Also, in order to make the euclidean Einstein action bounded from below, a "stabilized" 
euclidean action has been recently proposed. We shall briefly consider the effects of these 
modified kinds of dynamics on the correlation functions in Section 5.4. 

5.3 Two-point functions. 

An important point to be clarified about the vacuum correlation functions in gravity 
is the following: "which are the most suitable field quantities to be correlated" ? We recall 
that the criteria for the existence of gravitational waves [Zakharov, 1973] are usually based 
on the propagation of the curvature, which is believed to be the most physical effect of 
gravitation. One is thus led to consider as first candidates the following quantities: 

- the Riemann tensor Rp* (x); 

- the Ricci tensor R fiv ,(x) = R% a „(x); 

- the curvature scalar R(x) = g^ u (x)R flu (x); 

- the "rotation matrix", or plaquette TZp(x) = R c p ilu (x)ij^ 1 ' ', where <j^ u is an infinitesimal 
surface around x (71 must not be confused with the "Mandelstam covariant" of Chapter 
2). 

The linearized parts of R, R^ v and Rp^ are respectively given by 
R L = d 2 h* - d a d^h a p- 

K» = \( d2 K» + d ^K - d„d a h a „ - d v d a h^)- 

Rapnv = i^Padphpv - dpdyhuv - d a d u hf3^ + d(id v h ap ). 
Let us consider the invariant correlations 
G R {D) = (R(x) R(x')) ; 
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G Ricci {D) = {R^{x)U^' v ' \x,x') R^{x')) ; 
G Riemann (D) = (R^ l/ (x)U^'»' v '(x,x') R%^, u ,{x')) Q - 

G Loop (D,a,a') = (n%x)U^,{x,x')n%{x')) Q . 

Here U is the matrix of the parallel transport introduced in Section 5.1, computed along 
the geodesic joining x to x'; D is the geodesic distance between x and x' . 

On a flat background and to order k 2 these functions are very simple, since the 
matrices U(x,x') reduce to identity matrices and the tensors retain only their linearized 
parts. It is easy to verify that for the first three correlation functions one has [M., 1992b] 

Gr ~ G Ricci ~ G Riemann ~ k 2 d 2 5 /L (x - x') + o(k 2 ). (5.12) 

Thus these correlation functions vanish if x and x' are distinct, and they are not indicative 
of the correlation lenght of the system - at least in this regime of weak fields. 

The perturbative result above can be related to a general property of the correlation 
functions [see also Collins, 1984]. Let us consider a functional integral of the form 

z = J d[(f)(x)] e~ s ^ x ^ (5.13) 

and perform an infinitesimal arbitrary translation 8<p(x) of the integration variable <p(x). 
Since the integration measure is invariant with respect to such a transformation, we have 



J d[cf>(x)} 



-S[<i>(x)] -J dx 5<t>(x) E(4>{x)) 



where E((f)(x)) represents the classical equations of motion. Expanding the second expo- 
nential and disregarding terms 0((p 2 ) we obtain 

, = /*(*)]e-^>> (!-/™™), 

from which, by comparison with (5.13), we have 

J d[<f>(x)] e~ s ^ x ^ J dxS</>(x)E(<f>(x)) = 0. (5.14) 

If we choose 5<j)(x) to be of the form e5 4l {x), we conclude that 

W(x)))o = 0; 
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performing in (5.14) one more translation of the form e5 4 (x') we obtain by similar argu- 
ments 



(E(cf>(x)) E{<j>{x'))) = 0, 

and so on. That is, in a quantum field theory the vacuum correlations of field equations 
vanish. 

In the case of Einstein gravity, the equations of motion are 

Rpu(x) = 0; 

thus we have to any order 

(R^(x) R pa (x')) = 0. (5.15) 

Since in our perturbative approximation for the functions G the indices are contracted, in 
practice, just by 5^ u instead of the full metric g^, from (5.15) follows (5.12). 

Going back to the first-order analysis, we notice that the last correlation function, 
Gloo P , is identical, in this approximation, to a Wilson loop computed along a dumbbell- 
like contour (fig. 2). The evaluation of the Wilson loop will be the subject of the next 
Section, where we shall see, however, that it vanishes in this approximation. In conclusion, 
the only gravitational correlation functions which do not vanish to lowest order are those 
of the metric or of the vierbein 

G M etHc(D) = {g^{x)W^' v \x,x') gtx , u ,{x')) Q ; 

Gvierbein(D) = (e% (x) (x, x') e% (x')>0, 

as can be easily checked by a computation similar to that of the correlation functions of 
the curvature. They both behave like (compare also Section 4.4) 



5.4 Wilson loops. 



In the usual gauge theories the Wilson loop of the connection is one of the most 
important observable quantities ([Kogut, 1979]; see also the Section on the static potential 
energy, 6.4). 
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Unfortunately, in quantum gravity the Wilson loop of the Christoffel connection van- 
ishes, to lowest order, along any contourf . The proof of this surprising property is straight- 
forward: consider eq.s (5.7), (5.8) and denote by roman letters the corresponding vacuum 
averages. For instance: 

f/=(W)o = H(W (1) )o + ^A + ... 

= l + [ /(D + ^(2) + ... 

W = (W) = -4 + TrU 

= TrU (1) + -TrU {2) + ... 

= w w + \w^ + ... 

The Wilson loop W is given by 

W = \ P jdx» j dy v (T^(x)T^M) + o(k 2 ), (5.16) 
where the brackets denote the bare propagator of V. Using (5.2), (5.11) we obtain 

= -jk~ 2 { ai ^ r ^^ 1 ^ + aA '-' 64{x ~ v) }- (6 ' 17) 

where a\ and ci2 are two numerical coefficients. From (5.16), (5.17) one concludes that 
W vanishes to this order, except for some divergent perimeter contributions due to the 
5-function. 

In the next chapter we shall illustrate in detail the physical meaning of this vanishing 
and its consequences. About the causes, it is essentially due to symmetry reasons and to 
the absence of massive propagating modes. If one writes down the most general form a 

f Hamber ([1993]; see Chapter 7) argues that also in the non-perturbative sector the 
Wilson loop "does not have the same interpretation as in gauge theories, since it is not 
associated with the newtonian potential energy of two static bodies. In ordinary gauge 
theories at strong coupling the Wilson loop decays like the area of the loop, due to the 
strong independent fluctuations of the gauge fields at different points in spacetime and 
ensuing cancellations. In lattice gravity the situation is quite different since the connections 
cannot be considered as independent variables, and the fluctuations in the deficit angles 
at different points in spacetime are strongly correlated." The reasoning we present here 
holds in perturbation theory and to lowest order. See, however, Section 6.4 for a physical 
justification of why a closed loop cannot give in gravity the same result it gives in the 
gauge theories with "charges" of opposite signs. 
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massless graviton propagator may have while respecting Poincare invariance^, and repeats 
the calculation of the Wilson loop, one still finds that it vanishes to order %. 

For instance, the situation in (R+R 2 )-gravity (see Section 7.1) is different. Here R 
and R^, unlike in Einstein gravity, can propagate (cfr. Section 5.3). There can be some 
no n- vanishing invariant correlations of the curvature and eq. (5.17) will be replaced by 

W = {"' {- (2,/ x -vr ■ 

where Dm is the propagator of a massive mode. Therefore, the term proportional to ci2 is 
not ultra-local any more. However, we recall that only the lightest (massless) states with 
spin two should dominate at large distances [see for ex. Hamber, 1992a]. 

As it is well known, Einstein's gravity written in the first order formalism is a gauge 
theory of the Lorentz group (i.e., the action is invariant under local Lorentz transforma- 
tions), but not of the whole Poincare group ISO(3, 1). A consistent gauge formulation can 
be obtained only introducing some auxiliary fields q a [Grignani and Nardelli, 1992]. 

So it is not possible to consider in (3+1) dimensions, like in (2+l)-gravity [Witten, 
1988], the holonomies of the Lie algebra valued connection 

A»{x) = el{x)P a + Tf{x)u ah , 

where P a and u a b are the generators of the translations and of the Lorentz transformations. 
The holonomies of may have more content than the holonomies of alone. For 
instance, it can be easily verified that the term 

Trjdx^j dy»(e;(x)P a e h v {y)P b ) = -25 ab j dx» j dy»(e°(x) e h v {y)) 

is not trivial to leading order, unlike the corresponding term containing the connection. 
However, this term does not respect the invariance of the action. 



| A non-vanishing contribution proportional to % may arise on a non-flat background 
[M., 1993 c]. 
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6. FURTHER PROPERTIES OF THE WILSON LOOPS. 



6.1 Geometrical and physical interpretation. 

As we have seen in Section 5.4, the Wilson loop of the Christoffel connection vanishes 
to leading order in Einstein's gravity. Also, we remind that this loop is equivalent to that 
of the gauge connection (see eq. (5.6) and the ensuing discussion). 

In order to understand the physical meaning of this vanishing, a group-theoretical 
analysis of the matrix U is needed. We shall see that in the euclidean theory the vanishing 
of its trace amounts to a very strong geometrical statement. 

Let us first consider, for illustrative purposes, the case of a Yang-Mills theory of the 
group 50(3). The gauge connection has the form 

A lx {x)=A i lx {x)L i - i = 1,2,3, 

where the matrices Lj constitute a representation of the Lie algebra of the group. In 
particular, to fix the ideas, let us choose the adjoint representation; in this case the matrices 
Li have elements (L^)™ (m, n = 1,2, 3), which are related to the structure constants £i mn 
of the group. The connection A^(x) performs the parallel transport of a 3-dimensional 
vector V n in the "internal" space according to the formula (compare eq. (5.1)) 

dV m = Afa) (L l )™V n dx fl . 

The vector is rotated during the transport, but its length remains unchanged. Let us 
consider the matrix 0{C) which describes the parallel transport along a closed curve C. 
O(C) is defined by a P-exponential, through a formula similar to eq. (5.3). Suppose that 
we take a vector V in a point P of C, and parallel-transport it along C, returning to P; 
let us denote by V' the new vector we obtain in this way. The vectors V and V have the 
same length, that is 

SmnV m V n = 5 rnn V' rn V' n 7 (6.1) 

but they differ by an angle 6>, which is related to the trace of 0(C). For small angles, we 
have, by a proper choice of the coordinate axes in the internal space 



0(C) = 




that means 

TrC(C) = 3-9 2 . 

More generally, we remind that the Lie algebra of SO (3) has just one Casimir invariant, 
namely the operator 

L = L-^ -\- L<2 ~~\~ L^. 
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This operator commutes with each of the L^s, so we can in general rotate our coordinate 
system as to have L 2 = L§, and the rotation matrix takes in this case the form (6.2), i.e. 
we have 

0(C) = l + 9L 3 + ^e 2 Ll + ... (6.3) 

Taking the trace of (6.3), remembering that TrL^ = and using the normalization condi- 
tion of the Lie generators 

Tr LiLj = —25ij, 

we find that 6 2 is the coefficient of the Casimir invariant in the expansion of the exponential. 

Next we come to consider the group SO (4). Intuitively, adding a new dimension we 
can make an independent rotation. Multiplying two 4-dimensional matrices similar to 
(6.2), the first representing a rotation by an angle 6i perpendicular to one plane and the 
second a rotation by another angle 9u perpendicular to another plane, we find that 

TrO(C) = 4-(6j + 6j I ) (6.4) 

Also we know that SO(4) = [SO(3)]i x [SO(3)]n and that we have two Casimirs now [see 
for instance Wybourne, 1974], corresponding to (L 2 + L 2 j), whose "eigenvalue" appears 
in (6.4), and (L 2 — L 2 j), which is not of interest in this case. 

The group SO (4) is the relevant one for euclidean quantum gravity. In fact, the geo- 
metrical interpretation of the matrix U(C) is the following. During the parallel transport 
of a vector V in spacetime, its length, given by 

\V\* = V a V b 8 ab = V»V V g liV {x), 

does not change. If we transport V along a closed curve C, returning to the starting point, 
we obtain another vector V' , which has the same length of V , and differs from it only in 
the orientation. Hence we have for any vector 

V a V b 5 ab = V' a V' h 5 ab =U a c {C)V c U h d {C)V d 5 ab , 
or, in matrix notation, 

U T {C)U{C) = 1. 

The matrix U belong then to SO(4) and its trace has the form (6.4). 

If the variance of the angles Oi and On is zero to order % (because vanishes), the 
angles themselves have to vanish identically in any configuration, that is 

U{C) = 1 for any C. 
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This is a very strong geometrical statement, as it implies that, still to order h, all the weak 
field configurations which effectively enter the functional integral 

z = J d[h] exp {-h- 1 S[h]} (6.5) 

have no curvature. In other words, the curved configurations - which possibly dominate 
in other regimes - are in this approximation totally suppressed. 

This unexpected situation should be compared with what happens, for instance, in 
a ordinary SO (3) or SO (4) gauge theory. In this case the leading term W^ 2 \C) does 
not vanish and the variance of the rotation angles is not zero to order K. For instance, 
if the curve C has the form of a rectangle of sides L and T, with L < T, the quantity 
— (^T) _1 log(6> 2 )o is the potential energy of two non-abelian charges kept at rest at a 
distance L each from the other. 

So the matrices of the parallel transport in the "internal" gauge manifold, considered 
configuration by configuration, are not equal to the identity matrix. Interpreting Ti as the 
temperature G of an equivalent statistical system, we see that when G grows from zero to 
some small value - such that we may disregard G 2 or higher orders - the Yang-Mills fields 
develop "localized excitations" , i.e. regions of various sizes where the Yang-Mills curvature 
is not vanishing. 

All this does not happen for the gravitational field, which remains essentially in a 
"flat" state. Such a picture also explains the absence in this approximation of any invariant 
correlation of the curvature (compare Section 5.3). 

The curved configurations to order Ti in the functional integral could have been in- 
terpreted in a natural way as "virtual gravitons", since they represent gauge-invariant 
excitations which are off-shell and localized in space and time. (Virtual "gravitons of the 
metric" do exist in perturbation theory, but they do not necessarily represent physical ob- 
jects.) Their absence is a very peculiar property of gravity, which also forces us to imagine 
a peculiar mechanism for the gravitational interaction (see Section 6.4). 

6.2 The dumbbell correlation function. 

In this section all the matrix elements of U will be computed to order k 2 , for a 
dumbbell-like contour (see fig. 2). To this order we have = 0, while is given by 

U% = p£dx»£ dy» (T^(x) TZp(y)), (6.6) 

where C is the dumbbell. is invariant to order k 2 , by virtue of (5.10). Although not 
invariant to higher orders, it constitutes an interesting example of two-point correlation 
function. 

Disregarding gradient terms in the integrand and ultra- local terms of the form 5 N (x — 
y), eq. (6.6) becomes, in the Feynman-De Witt gauge and in any dimension N > 2 

U% = \ 5 W S aa ' P£dx»£ dy v ([dyha^ix) - d a ,h^(x)\ [dph Y „(y) - d Y Kf3(y)]) 

1 2 3N — N 2 f f d d 1 , n . 
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It is immediate to verify that the trace of is an ultra-local term. Eq. (6.7) is remarkable 
under two aspects. First, it shows that the matrix is symmetric. Second, we see that 
it vanishes identically in dimension N = 3; this happens because the Riemann curvature 
in an empty space-time of dimension 3 is always zero (compare the discussion following 
eq. (5.12). 

When C is a dumbbell the integral breaks down in 16 parts (fig. 3), which may be 
denoted as follows 

t/ (2) (r, D) = u A a + ubb + u C c + udd + 

+ U A B + u B a + u A c + uca + u A d + u DA + 

+ U B C + UCB + U B D + U DB + U C D + U DC = 

= uaa{t) + u C c(r) + 

+ ubb (d) + u DD (d) + 2u B D(d) + 

+ 2u AC (r,D). 

In this equation we have put d = D — 2r and omitted the indices a and (3. The "transport" 
terms depending on d can be easily shown to vanish in the limit of zero-width 5 of the 
strip (this holds only for the linearized theory, which behaves like an abelian one). 

The most interesting case is that of r <^ D. We remind that when a vector is trans- 
ported around an infinitesimal surface dx^ A dx v ', its variation is given by the matrix 

MP = RP^dx»Adx", 

where is the Riemann curvature. Like in Section 5.3, let us denote by 1Z? the matrix 
divided by the area, in the case of the surface being a small geodesic circle. is a 
kind of "regularized curvature" ; the regularization is non-local and gauge invariant, at the 
linearized level. Thus the integral t/( 2 ^(r, D) may symbolically be written in this limit as 

U%(r, D) = 2ttV {(^(0)7^(0)) + (T^(0)H}(D))} . 

The first term on the r.h.s. is divergent but independent on D and the second one is the 
"dumbbell function" (fig. 4) 

(7^(0)7^00)) = lim '^^dn+P fj ds fj dt ^os(s-t)k N r (s,t,a), 

where = r/D and the non-zero matrix elements are 

0i(s,t,a) = 1 -Nk 2 (l -au) 2 ; 
(pj(s,t,a) = 1 - Nk 2 a 2 v 2 ; 
4>l(s, t, a) = <j)1(s, t, a) = —N k 2 av(l — au); 
0£(s,£,a) = l; a* = 3,4,... 
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with 



k = l/\/l-2au + 2a 2 w) 

u = cos s — cost; v = sins — sin t; w = 1 — cos(s — t). 

Expanding these elements in powers of a and integrating!, one finds that the first 
non-zero coefficients are those of a 2 . In this way, denoting by $ the integral of <fi with 
respect to s, t, one finally obtains 



<S>\ = ^N(l + N)(N-2); 

®1 = -\n{N-2)- 
$2 = $1 = 0; 

^ = - l -N\ ^ = 3,4,... (6.8) 

In eq. (6.8) one can verify one more time that the trace vanishes. Actually, this is the 
point where this fact was noticed for the first time. 

6.3 Geodesic corrections. 

This section will be concerned with a problem raised by the presence of a "dynamic 
metric" in quantum gravity. The idea of a dynamic geometry is inherent to General 
Relativity, but it could not be applied to the quantized theory as long as it was considered 
as a field theory on a fixed background. A recent realization of this idea can be found in 
the Montecarlo simulations of Regge calculus (see Chapter 7). 

Our attention still focuses on the loops of the connection. A natural question arises: 
how can they be defined in the absence of a flat background? 

Let us first consider the case of a classical gravitational field. We assume spacetime 
to be described by a manifold M, endowed with a metric g. Let us consider on M a 
closed smooth curve C; C must be defined in an intrinsic fashion, without reference to any 
coordinate system. In particular, it is important to specify the form and the size of C, as 
in the two following examples: 

(1) the circle of geodesic radius r (fig. 5). It is defined as follows: given a center point O 
and a plane through O, we go along the geodesic lines which start from O tangent to the 
plane and we stop when the invariant distance from O is r. The points we reach in this 
way belong to the circle. 

f This is a quite long calculation, which can be easily performed using for instance the 
symbolic program Mathematica. 
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(2) the (planar) "dumbbell" (fig. 6). It consists of two circles of geodesic radius r, whose 
centers are placed along a geodesic line, the second lying an invariant distance D apart 
from the first. The circles are joined by a strip of small width 5 (S — > 0). 

In a curved space with arbitrary metric the geodesic circles and the geodesic dumbbells 
may "look" quite different from those in the flat space. When we compute the vacuum av- 
erage of a quantity like W through a numerical simulation of Regge calculus, what happens, 
intuitively, is the following: the computer produces an arbitrary field configuration, with 
a probability weighted by the exponential of the euclidean action; on this configuration W 
is "measured", along a contour defined in a similar way as we did for the geodesic circle 
or the geodesic dumbbell; these steps are repeated for many configurations and finally the 
average of the results is computed. 

On the other hand, the attitude one takes in the perturbative calculations, like in 
Chapter 5, is different. One works on a flat background, looking at small fluctuations of 
the metric around the background, but disregarding the effect of these fluctuations on the 
definition of the contours. 

It is possible to reconcile the two views in the case of a weak field, by introducing 
"geodesic corrections" in the perturbative calculations on a flat background [M., 1993a]. 
We start recalling some useful properties of the normal coordinates Chosen an origin 

O, they can be defined in a neighbourhood of O imposing the condition 

err^(o = o. (6.9) 

As a consequence of (6.9), the lines of the form £ M (A) = Au M , where v is a fixed vector, 
represent geodesic lines through O. Furthermore, the invariant distance between a point 
£ and the origin is simply given by ^/^^d^. It follows that a geodesic circle of radius r 
and center O has the familiar form 

f(s) = ( r cos s,r sins, 0, ...); < s < 2tt. (6.10) 

Given any coordinate system {x^} with origin O, the normal coordinates £ can be 
computed from the coordinates x by solving the system (3.11), (3.12). We are interested 
in the first-order perturbative solution given by (3.13), namely 

x»(0 = e + x»(t) + o(r 2 ), (6.ii) 

where 

= -er f dt (i - 1) r^o, (6.i2) 

Jo 

being V the connection in the coordinates x. 

The corrections due to (6.11) involve in the term which on a flat background 

gives no contributions to order k 2 . We choose {x^} to be the harmonic coordinates fixed 
by the Feynman-De Witt gauge. In order to fix the ideas we also make the hypotesis that 
C is a geodesic circle, although this hypotesis is in fact irrelevant. The function x^(s) 
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which represents the circle in the coordinates x can be obtained from eq.s (6.10), (6.12). 
Since we are interested in terms up to order T 2 (corresponding to order n 2 in the quantum 
perturbative series), the approximation (6.12) will be sufficient. We thus have 

+ <fde>d v rfo(£)x"(t) +o(r 3 ), (6.13) 

where £ is given by (6.10). In eq. (6.13), the first integral represents the naive contribution, 
which vanishes when it is averaged on the vacuum; the second integral is a kind of "index" 
correction; the third one is a correction to the "argument". Using (6.12), one finds for the 
average value of U^ 1 ' on the vacuum state 

U% = - f dt(l-t)^{t)+o(K% (6.14) 
Jo 

where 

fg(t)= I d(£'n(W0rW<0>+ 

J c 

+t £dc ^(r^(Od^ pa (to)+ 

Disregarding gradient terms in the integrand we obtain 

/|(f) = 6°"' ~S% I df £T (PA W0] [o 'fl T^M)) ■ (6-15) 

The antisymmetry of fp causes its trace to vanish. It follows that the geodesic contri- 
bution of order k 2 to W vanishes, like the naive contribution TrU^- 2 \ We notice that the 
antisymmetry of / in (6.15) does not depend on the form of the propagator, and is thus a 
geometric property. On the contrary, we remind that the trace of the naive contribution 
is sensitive to the dynamics of the field (compare Section 5.4). We finally point out that 
eq. (6.15) holds for any closed curve C, not just for the geodesic circle. 

Returning to (6.15), let us now compute all the matrix elements. We notice that eq. 
(6.15), substituted into (6.14), produces a correction of the form 

U {1) = f dt (1 - t) <f d£ (*(£) *(*£)), (6.16) 
Jo Jc 
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where \I/ is a field. This expression shows that the geodesic contribution is related to the 
necessity of fixing a scale for the distances. Since (6.16) diverges for t ~ 1, in the following 
we shall regularize it by stopping the integrations at t = 1 — e. 

Substituting the propagator (5.11) into (6.15) we find, after some steps, in any dimen- 
sion N > 2 



U% = -c N £ dC *T x 



[a'0\, (pa) 



N — 2 J ' ' \pf_2j tXT p f,a £N-2 

Finally, disregarding gadient terms in the integrand we obtain 

U% = ~c N k 2 X #(e) £ \ [dCdp - d^d a ] = 

= -cn xP (e) £ d(C A fa -L (6.17) 

where Xn\ 6 ) diverges like e 1_Ar when e — > 0. For a geodesic circle of radius r in the plane 
1-2, the only non zero elements of are simply 

For a geodesic dumbbell we have, in the limit r <^ D, 

U^\ = -U^\^-c N ^ X %\e) 



2 

7rr 



(#)" 

This correction should be added to (6.7). Although it is formally quite clear, a further 
analysis seems to be necessary in order to understand its physical meaning. 

We also remind that the matrix (6.17) is not invariant under gauge transformations. 
Nevertheless, its simplicity suggests us that it could have a quite general geometrical 
meaning. 



6.4 The static potential energy. 

It is known that the static potential energy U(L) of two sources of a gauge field is 
related to the Wilson loop of temporal size T and spatial size L « T by the formula 

e -h-*TU{L) =W (L 1 T). 

In the strong coupling limit this expression leads to the confining potential U(L) = kL, 
while in the weak coupling limit one easily recovers the Coulomb potential U(L) = —e 2 /L. 
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The vanishing of W to leading order in gravity raises the problem of finding another 
invariant expectation value of the quantized field which gives the static potential energy 
between two masses. This problem has a well defined solution indeed [M., 1993b]. One 
starts from the following known formula of Euclidean field theory 

£ = lim i fmeM-n-\sM + s lnte ^j ])} 

t-oo %T s Jd[(f>] expf-n-^oM} 

where is a quantum field and J is a classical source coupled to 0, which is switched off 
outside the interval (—^T, \T). £ represents the energy of the ground state of the system. 
We shall only show here how this formula works in the case of a weak gravitational field 
on a flat background. Replacing 4> with the gravitational field, So with Einstein's action 
and J with two particles of masses mi, m 2 , following the trajectories 

= ~, 0, ; y»(t 2 ) = (t 2 , |, 0, , (6.19) 

we find (for a weak field) 

£ = lim -J=x 

T^oo hT 

J d[h] exp h~ x I -SWnst. [h] - mi J dti a/1 - h 00 [x(ti)] - m 2 J dt 2 a/1 - h 00 [y(t 2 )] \ 

x log -. 

J d[h] exp{-^ SEiwt.[h]} 

It is easy to verify that to leading perturbative order this gives the correct result 

m x m 2 G 

£ = mi + m 2 j • 

It is interesting to make a comparison with electrodynamics. In that case the analogue 
of the functional integral which appears in the logarithm of (6.18) has the form 



2 C 2 



2 2 



exp <^ e / dhAoixih)] - e / dt 2 A [y(t 2 )] } . (6.20) 




(The two charges have been chosen to be opposite: q± = e, q 2 = — e.) Reversing the 
direction of integration in the second integral and closing the contour at infinity, one is 
able to show that the quantity (6.21) coincides with the Wilson loop of a single charge g, 
thus giving a gauge invariant expression for the potential energy. 

In gravity this is not possible: we may imagine that an expression like (6.21) could be 
obtained in the first-order formalism (with A replaced by the tetrad e[j), but the masses 
necessarily have the same sign, so the loop cannot be closed. 

It is possible to give a definite meaning to the expression (6.18) for the energy also 
in the case when no background is present. It can also be proven that the gravitational 
interaction energy is always negative [M., 1993b]. A very recent work based on this formula 
is that of Hamber and Williams [1994]. 
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7. LATTICE GRAVITY. 



For quantum field theories with relevant non-perturbative effects, the lattice discretiza- 
tion is often quite important, since it allows numerical evaluations through the Montecarlo 
simulation technique. This is also the case for quantum gravity, which is usually thought 
to be very strong at small distances, where the theory should somehow determine the small 
scale "quantized" structure of spacetime [see Rovelli, 1991 and ref.; Smolin, 1992 and ref.]. 

In the discretized theory the correlation length is a crucial quantity. For this reason 
we shall briefly recall here some advances made in lattice gravity during the last years, 
especially from Hamber, Caracciolo and Pelissetto, Ambj0rn and co-workers and Greensite. 

The point of view of the lattice [see for instance Kogut and Wilson, 1974; Kogut, 
1979; Seiler, 1982; Creutz, 1983] is very different from that of continuum field theory; the 
two are connected, however, by the "continuum limit" . 

Let us consider the equivalent statistical system of a field theory discretized on a 
euclidean lattice. We remind that if the transition between two phases of the system is of 
the second order, then the appearence of long range correlations causes the details of the 
lattice to become irrelevant. The original (renormalizable) continuum theory can thus be 
reobtained, at the transition point, when the lattice spacing goes to zero. If the original 
continuum theory is not renormalizable in the usual sense (like Einstein gravity), a more 
general procedure, called "asymptotic safety" [Weinberg, 1979] can be possibly applied. 

The chapter is organized as follows. As schematically shown in Table 1, the study of 
a field theory on the lattice starts choosing one specific version of the theory and a dis- 
cretization scheme. Some possibilities for gravity are displayed in the table for illustration 
purposes, with no attempt to completeness; in fact, we are not interested here to make 
a review of all the different approaches and to compare the results (which should be the 
same, at least for physical quantities!); we are interested here to the vacuum correlations at 
geodesic distance. So we have chosen to analyze in more details the way which seems to be 
the most promising under this respect, namely: (R + R 2 ) -gravity (Section 7.1), discretized 
through Regge calculus (Section 7.2), as it appears in the simulations of Hamber (Section 
7.3). The remaining section, 7.4, is dedicated to a brief description of other approaches. 

Throughout the chapter we shall often quote from the review of Menotti [1990], were 
many omitted details can be found. 



7.1 (R + R 2 )-gravity. 



It is known that the scalar curvature R is the only scalar which contains second 
derivatives of the metric. Further scalars can be constructed, however, which contain higher 
derivatives, namely R 2 , R^ U R^ V , R^ upa ' R^ vpa (among the three, only two are independent). 
The Einstein's action can thus be modified as follows: 



S 



-I 



1 R-^R^ upa R^ 



K 



(7.1) 



The terms quadratic in the curvature would not be observable in macroscopic phenomena, 
but they are important to the quantum theory, since the first radiative corrections to 
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Einstein's lagrangian contain just R 2 terms. As Stelle showed in detail [1977], (R + R 2 )- 
gravity can be considered as a renormalizable, although non unitary, theory f. Various 
attempts of recovering unitarity by graphs resummation, suitable gauge fixing and other 
techniques [Antoniadis et al., 1986; see also Weinberg, 1979] did not lead to any definitive 
outcome. 

The authors of numerical simulations on (R + R 2 ) models believe usually that the loss 
of unitarity is a problem confined to perturbation theory, as the whole theory could still 
be unitary [Hamber, 1986]. 

7.1 Regge calculus. 

The basic idea of this formalism [Regge, 1961] is to replace the smooth spacetime 
manifold with a simplicial manifold, i.e., a manifold made up of flat regions ( "simplexes" ) 
connected by "edges". This is easy to figure out in two dimensions (see fig. 7). The 
discontinuity in the parallel transport arises when an edge is crossed, and the curvature 
concentrates in those points where a few edges meet, called "hinges". If we transport a 
vector along a closed path, we only notice a difference between the original vector and the 
transported one, if the path contains at least one hinge. In a generic dimension A the 
dimension of the hinges is A — 2; hence for A = 2 they are points, but for A = 3 they are 
lines and for A = 4 they are triangles. 

Given an hinge and one simplex which touches it, the opening angle On of the simplex 
with respect to the hinge can be defined. In two dimensions, the geometrical meaning of 
such an angle is apparent. In general, n is given by the expression 



. A V N V N -2 

sin pat = 

A — 1 Vn-i V'n-i 



where Vn is the volume of the mentioned A-dimensional simplex and V/v-i, V'n-i are 
the volumes of the two (A — l)-dimensional simplexes which share the hinge (in the two- 
dimensional case, they are the two edges which meet at the hinge). 

Let us now consider all the simplexes which meet at a given hinge and sum their 
deficit angles referring to that hinge. We obtain in this way a "deficit angle" 

sCh 

Regge postulated the following, simple action: 

S = 2Y,A h 6 h , (7.2) 

h 

where the sum ranges over all hinges and A h is the surface of the hinge. In the contin- 
uum limit, the action (iliz) really approaches Einstein's action (for references about the 



f It has been shown that in a suitable range of the couplings a and b the theory is 
asymptotically free with respect to them. Also, asymptotic freedom with respect to G 
may arise under certain circumstances [Julve and Tonin, 1978] 
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exact meaning of this limit and the uniqueness of Regge's action, see [Menotti, 1990] and 
[Williams and Tuckey, 1991]). 

The extension of Regge's action by a cosmological term of the form A Yl s V s is straight- 
forward. On the contrary, the lattice version of the R 2 terms is quite involved [for a review 
see Hamber, 1986]. To write them, one must associate to each hinge a suitably defined 
volume Vh- The simplest R 2 term is 



The other terms are more complex, as they entail the introduction of a lattice analogue 
of the Riemann tensor. All the lattice translations are not unique; however, they are 
equivalent in the continuum limit - when the maximum lenght of the edges vanishes, with 
respect to the characteristic radius of curvature of the smooth manifold. 

7.3 The simulations of Hamber. 

In a remarkable series of papers starting in 1985, Hamber has developed a technique 
for simulating quantum Regge calculus of (R + i? 2 )-gravity [see Hamber, 1986, 1992, 1993 
and references therein]. 

The typical quantities which were first computed in these simulations were rather 
reminiscent of Ising model, namely the average curvature R and the average curvature 
squared R 2 . It turned out that they diverge along a line which goes across the parameters 
space (a, «), thus dividing it into two phases, called "smooth phase" and "rough phase". 
The critical indices were computed too. 

Recently, the evaluation of two-point correlation functions at geodesic distance has 
been addressed [Hamber, 1993 b]; the results are interesting although rather preliminary. 

In a typical simulation, the simplicial lattice is built up of rigid hypercubes, which can 
be subdivided into simplices by introducing face diagonals, body diagonals and hyperbody 
diagonals. This choice is not unique and is dictated by a criterion of simplicity, with 
the advantage that such a lattice can be used to study rather large systems with little 
modification. The length of the edges / is individually varied (by moving at random 
through the lattice), and a new trial length is accepted with probability min(l, e _AS ), 
where AS is the variation of the action under the change in edge length. If the triangle 
inequalities or their higher-dimensional analogues are violated, the new edge length is 
rejected. It should be noticed that in order to compute the variation in the action under 
the change of one edge length, a large number of adjoining triangles and their deficit angles 
has to be considered. 

Lattices of size between 4 4 (with 3840 edges) and 16 4 (with 983040 edges) were con- 
sidered. Periodic boundary conditions (four-torus) were used, since it is expected that for 
this choice boundary effects should be minimized. The lengths of the runs typically varied 
between 10-40/c Montecarlo iterations on the 4 4 lattice, 2-6/c on the 8 4 lattice and 0.5k on 
the 16 4 lattice. On the larger lattices duplicated copies of the smaller lattices are used as 
starting configurations, allowing for additional equilibration sweeps after duplicating the 
lattice in all four directions. 
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Like we said above, typical quantities which are "measured" in the simulations are 
the average curvature 



2x (2J2 h 5 h A h ) (J^gR) 



R {l/ (Eh V h) " (IV9) 
and the average of the curvature squared 



(Eh Vh) (J V9) 

One can also estimate the local fluctuation 



Xr 



(EhVh) 



U^hAh) >-<2 5>A) 

V h J h 



A divergence in x should be indicative of a second-order phase transition (as it happens 
for the magnetic susceptibility or the specific heat). 

As the bare Newton constant and the coupling a are varied (compare eq. (7.1), where 
one sets in this case b = a/4), a continuus phase transition is found, separating a "smooth" 
phase from a "rough" phase. In the first phase the curvature is small and negative, and the 
fractal dimension is consistent with four. In the second phase the simplices are collapsed, 
the curvature is large and positive, and the fractal dimension is much less than four, 
indicating the presence of finger-like structures. Approaching the critical point from the 
only physically acceptable phase, the smooth one, it was found that the curvature vanishes 
with an exponent 5 = 0.62(5). At the critical point the curvature fluctuation x diverges, 
leading to the possibility of defining a non-trivial lattice continuum (this can be then 
checked measuring the correlations: see below). 

We notice thus that the vacuum expectation value of the curvature can be used as 
an order parameter for the transition; moreover, it can be used as a possible definition of 
the effective, long-distance cosmological constant. Usually one adds to the action (7.1) a 
cosmological term with some bare A; the effective cosmological constant is then given by 

4A\ _ (fygR) 

K Jeff (IVd) 

As one approaches the fixed point at k c (with fixed a), one finds as expected (A/ft) e ff — > 0. 

The computation of the correlations at geodesic distance has been recently addressed, 
in particular for the averages 

G R (D) ~ (J2 S h A h J2 5 h' A h> S(\x - y\ - £>)), 

h(x) h'(y) 

which corresponds to the correlation of the scalar curvature, namely 

~ (y/0(x)R{x) Vg¥)R(y) s(\x -y\- d)) 
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and 



G V (D) ~ E v * 6 (\ x -y\- D ))> 

h(x) h'(y) 

which corresponds to the volume correlations 

~ (VfKx) y/g(y) s (\ x ~v\ ~ D ))- 

The computation proceeds as follows. First, the geodesic distance between any two 
points x and y is determined in a fixed background geometry. Next the correlations are 
computed for all pairs of points within geodesic distance D and D + AD, where AD is 
an interval slightly larger than the average lattice spacing Iq = a/ (l 2 ), but much smaller 
than the distance D considered. Finally, the correlations determined for a fixed geodesic 
distance D are averaged over all the metric configurations. 

In principle one could also compute correlations of vector and tensor quantities in- 
troducing the matrix of the parallel transport (compare Chapter 5), but this is quite 
complicated on a Regge lattice and has not be done yet. 

The details of the method are described in the mentioned paper. The main result is 
that the volume correlations are negative at large distances, while the curvature correla- 
tions are always positive. If the correlations are fitted to an exponential decay, one finds 
that the behavior is always consistent with a mass that decreases when one approaches 
the critical point. This behavior is precisely what is expected if the model is supposed 
to reproduce the classical Einstein theory for distances which are very large compared to 
the ultraviolet cutoff scale. Note that this happens for a model of gravity which at short 
distances is far removed from the pure Einstein theory, containing both a bare cosmological 
term and bare higher derivative lattice terms! 

7.4 The gauge approach and others. 

We shall describe the main features of the gauge formalism in dimension N = 4, 
referring to the literature for issues like the reflection positivity, the invariance under 
reparametrization transformations and the phenomenon of graviton doubling (see Menotti 
[1990] and references; Smolin [1979]; Menotti and Pelissetto [1986, 1987]). 

In the gauge approach one exploits the (incomplete) analogy between gravity and the 
gauge theories. The attitude is taken of considering gravity as a field theory in Minkowski 
or euclidean space, whose lagrangian is invariant under a group of local transformations 
(reparametrization transformations). Several people have introduced along these lines 
discretized versions of gravity. We shall briefly describe here the formulation due to Smolin, 
which can be considered as the prototype of these gauge formulations. Essentially Smolin's 
formulation is the discretized version of the Mac Dowell - Mansouri gauge formulation of 
De Sitter gravity [1977]. Let us consider the De Sitter group 0(4, 1), which goes to 0(5) 
in the euclidean, and introduce the usual hypercubic lattice, familiar to gauge theories. 
Associate now to each link of the hypercubic lattice a finite element of 0(5) 

U(n,n + n) = exp (|j bc wj c + aP c e c ^j , 
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where J a b and P a satisfy the following commutation relations 




SacJbd — SbdJac + SadJbc 



The action takes the form 



where a = 



3 



and 



Up?(n) = U AC {n,n + y) U CD (n + //, n + fi + v) x 
U DE (n + n + v,n + v) U EB (n + v, n) 



is the plaquette in the vector representation. It is clearly invariant under local 0(4)- 
transformations, but not under 0(5). This is a consequence of the incomplete similarity 
between gravity and gauge theories. In the formal continuum limit we obtain the Mac 
Dowell-Mansouri action 



which is the familiar Einstein action, plus a cosmological and a Gauss-Bonnet term. The 
discretized action is reparametrization invariant only in the formal continuum limit. 

The extension of the procedure above to the Poincare group has been given by Menotti 
and Pelissetto [1987]. To this end one has to introduce in the formalism the local tetrads, 
which transform locally under 0(4). They are defined by the expression 



where H = diag(l, 1, 1, 1, —1). 

The gauge approach is strictly bounded to the program outlined by Weinberg of finding 
a fixed point around which the theory could possibly be renormalized. 

The simulations were performed by Caracciolo and Pelissetto [1987, 1988; see also 
Menotti, 1990] on the De Sitter model in the gauge formulation. They made extensive 
simulations on an 8 4 lattice. Due to the compact nature of the group SO (5) and the 
presence of the lattice cut-off, the action is limited and thus the theory must exibit a well 
defined ground state. For small values of a the vacuum is similar to the QCD vacuum 
in the strong coupling region. The large a region is characterized by the vanishing of the 
vierbein, so that the action becomes dominated by the topological term. An interesting 
order parameter is given by 




E a (n) = --J-tr [HP a U(n,n + fi)HU(n + n,n)\ , 



P = 
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whose distribution is peaked around zero in the small a region and around the extremes 
±1 in the large a region. At the value ^ = 0.08 ± 0.01 a sudden phase transition is 
observed. Around this value all measured quantities, like the mean value of the action 
per site, the trace of the metric, the 0(4) curvature and P, exibit very large hysteresis 
cycles, thus showing that the transition is first order. For a > a c the system ends up 
in non reproducible states showing the existence of many metastable states above a c and 
recalling a situation similar to a spin glass. 

These results are very different, also qualitatively, from those of Hamber. Unfortu- 
nately, such disagreements are not unusual in the numerical simulations of very complex 
models. In a third approach [Ambj0rn et al., 1992] a transition is found, that could be 
first- or second-order; furthermore, the average curvature at the transition is non zero, 
which makes the continuum limit quite difficult to understand. 

Finally, it is worth mentioning that recently Greensite has given a new interesting 
approach to euclidean quantum gravity [Greensite, 1993 and references]. It is known 
that the euclidean Einstein action is not bounded from below, and this means that the 
euclidean functional integral is not well defined. There have been various attempts to 
solve this problem [see for instance Mazur and Mottola, 1990 and references]. The method 
of Greensite is in practice a prescription which can be related to stochastic quantization 
and which prevents the field configurations in the functional integral from running into 
the singular configurations which make the action unbounded. This method leads to an 
effective action which reproduces Einstein's equations but is non-local beyond first order. 
Although higher order calculations are difficult, it can be proved in a nice way [Greensite, 
1992] that the effective cosmological constant generated by this model is vanishing. The 
Montecarlo simulations are still at an initial stage. 
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